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ABSTRACT 


In this report electrical methods of determining soil 
moisture content are explored. Since the magnetic permeability 
and electrical conductivity of soils are known to be unreliable 
indicators of soil moisture content; the report focuses on the 
electrical permittivity of soils. The first part of the report 
gives an assessment of permittivity as an indicator of soil 
moisture content, based on experimental studies performed by the 
authors. The conclusion is that the electrical permittivity of 
soils is a useful indicator of aviilable soil moisture content. 

In the second part of the report, two methods of determining 
the permittivity profile in soils are examined in light of the 
findings in Part I of this report. A method due to Becher is 
found to be inapplicable to this situation. A method of Slichter, 
however, appears to be feasible. The results of ^lichter's 
method are extended to the proposal of an instrument design that 
could measure available soil moisture profile (percent available 
soil moisture as a function of depth) from a surface measurement 
to an expected resolution of 10 to 20 Cm. Extension of the results 
to the airborne , remote sensing problem is considered. 
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INTRODUCTION 

In the field of agronomy, many scientific investi- 
gations, as well as field applications, make use of soil 
moisture content data. Soil moisture is also an important 
variable in water resources management. Unfortunately, 
however, soil moisture content is not easy to measure in 
the field. The commonly used methods (gravimetric, con- 
ductivity, neutron thermalization, gamma ray attenuation) 
all have important disadvantages. Thus, there is a need 
for a better technique which can be used in remote sensing 
applications to field work. The method should have better 
accuracy than the conductivity technique, greater speed 
than the gravimetric method, and better portability and 
safety than the radiation techniques. 

Ideally, the method should permit the determination of 
the available (for plant growth) soil moisture profile (i.e., 
% available moisture as a function of depth). Electrical 
methods of measuring soil moisture appear to be able to sat- 
isfy these needs (32). Three properties of matter can be 
measured by electrical methods: magnetic permecibility , con- 

ductivity, and electrical permittivity. Except under un- 
usual circumstances, the permeability of soil is very close 


to its value in free space, Uj (1). The conductivity of 
soil is not a reliable indicator of moisture content be- 
cause of the large influence of ions in the soil. Mea- 
suring the conductivity of porous blocks buried in the 
soil overcomes this difficulty, but leads to inconvenience 
and inaccuracy. Only permittivity remains. It is plaus- 
ible that permittivity could be a reliable indication of 
soil moisture, since the dielectric constant (relative per- 
mittivity) of soil is around five while water has a relative 
permittivity of eighty. 

Furthermore, the electrical permittivity method of 
soil moisture measurement discussed herein has three im- 
portant advantages. First, it is well suited to remote 
sensing, which is important from the standpoints of economy 
of time and financial resources required. Second, it has the 
capability of measuring the moisture content as a function 
of depth below the surface of the earth, to a depth of the 
order of a meter. Third, it is not affected by irregularities 
in the flatness of the soil, of the order of a few centimeters 
within a radius of a meter, or so. 

It is possible, of course, to determine the amount of 
moisture at various depths in the soil (i.e. , the moisture 

I 

profile) by augering a hole and lowering some sort of probe 
into the soil. Another technique is to remove a core of the 
soil and to measure the soil moisture directly. Both of these 
methods are very time consuming. The proposed electrical 
permittivity method would not disturb the soil in any way since 




3 


the electromagnetic waves, to be used (at a frequency of about 
300 megahertz (MHz) or a free-space wavelength of one meter) 
can penetrate most soils to a depth of about a meter and al- 
so determine the available soil moisture profile with a def- 
inition in soil depth of the order of ten centimeters, or so. 

Active cutd passive radar methods are under consideration 
by other orgemi rations. These normally involve airborne in- 
struments operating at frequencies ranging from a few giga- 
hertz (GHz) to perhaps 25 GHz (33). In general, these methods 
(operating at wavelengths from about 1 to 10 centimeters) are 
handicapped by the existence of earth surface roughnesses of 
the order of a few millimeters and by the fact that, at these 
high frequencies, the waves penetrate the earth no more than a 
few centimeters, at the most, depending upon the conditions 
(moisture content and free ion content) of the soil. Further- 
more, they give a measure of the total moisture content of the 
soil, from the surface to their deepest point of penetration. 
This is a serious limitation for several of tlie more important 
applications of the measurements. The use of lower fre- 
quencies (longer wavelengths) is limited by the required in- 
crease in size of the antennas to be carried by the aircraft. 

The research involved in this present project naturally 
divided itself into two parts. The first part was concerned 
with determining accurate quantitative relations between soil 
moisture and electrical permittivity for various common soils, 
and the second part dealt with methods for determining the per 
mittiyity profile. The first part of this research has been 
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completed and the results are contained herein. The inves- 
tigation of the second part seems to indicate that the most 
promising method for meeting the end requirements involves 
the development of em electrical instrument to be placed on 
the surface of the earth (31). The measurements made by this 
instrument on the soil moisture profile could then be tele- 
metered to some central point, if desired. As mentioned pre- 
viously, it. is anticipated that the electrical measurements 
would be made at a frequency of about 300 MHz (a free-space 
wavelength of cLbput one meter) . 

Since completing the portion of the project described 
in the present report, cm additional report (3^) has been re- 
ceived. It discusses many methods for measuring the permit- 
tivity of the earth, but none seems to be particularly ap- 
plicable to the present need. It does contain, however, an 
excellent list of references on the subject. 
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PART I 

The Relationship Between Soil 
Moisture and Electrical Permittivity 
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CHAPTER 1 

STATEMENT OF THE PROBLEM 

In part one of this thesis, the relationship between the 
electrical permittivity of soils and the moisture content of 
soils is explored. 

Objectives of the Project 

We wish to determine whether permittivity is a good indi- 
cator of soil moisture content. Certainly, moisture content in- 
fluences soil permittivity. The fact that soil permittivity is 
a function of moisture content, however, is not sufficient to 
conclude that permittivity is a good indicator of soil moisture. 
It is necessary to explore other factors which affect soil per- 
mittivity and to assess their roles. 

The objective of the study is to answer the question; Is 
it possible to determine soil moisture content with acceptable 
accuracy from knowledge of the permittivity alone? 

Definitions 

Since this report is an interdisciplinary investigation, it 
may be helpful to all to clarify some unfamiliar terminology. 

Soil Moisture Content is the amount of water in the soil, 
generally expressed as a percentage. 
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1-1 % Moisture content = Mass of Water in Soil 

Mass of Dry Soil 

Wilting Point Moisture is the soil moisture content at which 
plants undergo permanent wilting. Practically speaking, it is 
defined as the percentage moisture content corresponding to a 
soil moisture tension of 15 bars. 

Field Capacity is the soil moisture content of soil which 
has been thoroughly wetted, but all the moisture that can drain 
off due to gravity has done so. This generally corresponds to a 
soil moisture tension of one-third bar. 

Available Water is the amount of water between the wilting 
point and field capacity. For a given moisture content, the 
percentage of available moisture is (Percentage Moisture - Per- 
centage Moisture at Wilting Point) . . 

Soil Moisture Tension is the force per unit area which binds 
the water to the soil. This is measured by the "suction" 
(pressure difference) required to extract moisture from the soil. 
The usual technique for obtaining a certain moisture tension is 
to place a soil sample next to a porous membrane in a pressure 
vessel and apply air pressure to the inside of the vessel. The 
pressure differential removes any moisture held by the soil par- 
ticles at less tension than the applied pressure. 

Permittivity (e) is an electrical property Of materials de- 
fined by Coulomb's Law. Rationalized MKS, units are used for all 
physical quantities. 

1-2 F=qiq2/4TT er^ 

where F is the force between two charges. 
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qj and q 2 are the magnitudes of the two charges, 
r is the distance between the two charges, and 
e is the permittivity.’ 

E can be expressed as where is the permittivity of a 

vacuum in C^/N-m^. e^ is referred to as the relative per- 
mittivity, or dielectric constant. 

Dielectric Constant is defined to be the relative permitti 

vity. 

Permeability (yi) is the ratio of magnetic induction to 
magnetic field strength (for linear materials) . In vacuum, the 

~ 7 

permeability is defined to be 4tt x 10 henry/meter. In this 
thesis, it is assumed that the permeability is equal to the 
vacuum value. This is almost exactly true, except for ferrous 
materials, which do not occur in the present work. 

Conductivity (a) is the ratio of volume current density to 
electric field strength. Conductivity can be calculated quite 
easily from resistance measurements. 

Some Factors Which May Affect the Permittivity of Soil 
Many factors influence the permittivity of soil. Some of 
the most important ones are introduced below. 

Soil Types. Soils are composed of many different types of 
materials whose dielectric cOhstants may differ considerably. 
Furthermore, soil and water interact. It is by no means incon- 
ceivable that different types of soil may interact differently, 
so that the permittivity of the mixture may differ significantly 
among different soils. 

Soil Condition. It is possible that the physical condition 
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of the soil may affect its dielectric properties. Compaction is 
probably the most important factor. A rather highly compacted 
soil has less airspace than a less compacted one. This might 
have some fund^ental effect on the permittivity of the mixture. 

In addition, compaction may have some other consequences. For 
example, suppose permittivity were proportional to the number of 
water molecules within a sample holder. If two samples had the 
same percentage moisture content, a more compacted sample would 
have more molecules of water and, therefore, have a higher die- 
lectric constant. 

Ion Concentration. The ion concentration in soils has a 
very strong effect on the electrical conductivity of soils, but 
it is not apparent that ion concentration will influence permitti- 
vity. If ion concentration is an important factor, this fact 
could have serious consequences. Ion concentration varies from 
place to place and time to time. Irrigation practices and fer- 
tilization alter ion concentration. 

Time. A number of references have been found in the litera- 
ture to the variation of measured permittivity with time. Var- 
ious materials were involved, including vanadium oxide. No ex- 
planations were given except that the investigators felt that 
these observations were not artifacts. 

Temperature. Temperature affects the permittivity of most 
substances. It is known that the permittivity df water decreases 
with increasing temperature. 

Frequency. The permittivity of any material is a function of 
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the frequency of the electromagnetic fields used in the measure- 
ment. A material- whose permittivity does not change over a 
certain range of frequency is said to be nondispersive in that 
range. 



CHAPTER 2 


LITERATURE REVIEW 

For several decades dielectric methods of determining 
soil moisture content have attracted consider£03le attention. 
The first section of this chapter will review some of the 
papers in this area, while the second section will deal with 
specific aspects of soil moisture interaction. 

Determination of Soil Moisture 
Using Dielectric Methods 

The earliest article to be surveyed is "A Review of Re- 
sults of Dielectric Methods for Measuring Moisture Present in 
Materials” by N.E. Edlefsen (2) . He discussed a number of 
moisture meters, mainly for food products. However, he did 
discuss three investigations devoted to soil moisture measure- 
ment, the investigators cited being W.L. Balls (3,4), G.H. 
Cashen (5), and the author himself. Balls reported a linear 
relationship between capacitcmce and moisture for all mois- 
ture contents except very low ones, using buried electrodes. 
Cashen used a different type of condenser and obtained ir- 
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regularly shaped curves which were difficult to explain in 
terms of soil moisture. Edlefsen also found a linear re- 
lationship between capacitance and water content, although 
he did mention difficulties in obtaining uniform compaction. 

He also reported successful field tests. Edlefsen continued 
his research, and in 1934 reported the development of "A Mew 
Capillary Potentiometer" (6) . 

During the 1930 's several investigators worked with di- 
electric properties, notably Anderson, Shaw and Muckenhirn, 
Aleksandrov, Banner jee and Joshi, and Yevstigneyev. Refer- 
ences to most of these can be found in Fletcher (7). in 1939, 
Fletcher published the results of his research. He built 
probes of plaster of Paris between metel electrodes, and his 
laboratory results were encouraging. Not only was he able to 
obtain a useful relationship between soil moisture and measured 
dielectric constant, but he was also able to shpv; that between 
0.1 and 20 grams of NaCl per 100 cc of water had no effect on 
his readings. 

The next study was that of Anderson and Edlefsen (8) who 
worked with an apparatus similar to Fletcher's. The invest- 
igation was designed to answer three basic questions: 

{1) Are the results obtained by dielectric methods re- 
producible? 

(2) How much time must be allowed for the block to come 
to moisture equilibrium? 

(3) What is the effect of plate separation? If electrode 
separation is critical, the expense of manufacturing 
moisture meters would be greater. 
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They found that the results were reproducible , the time 
lag was small if the probe was located near rapidly transpiring 
roots, and "no outstanding differences for electrode separations 
of 4 cm and 2 cm". 

Anderson gives the equation, C = which means that the 

capacitance of any capacitor is equal to its capacitance when 
filled with air multiplied by the dielectric constant of the 
material used to fill it. In the case of Anderson's sample 
holder, was no more than 0.0001 pF, the smallest value 
measured, and C = 0,1 yF, the largest. Then D = C/Cq = 1000. 
Soils have a dielectric constant of around five while water 
has a relative permittivity of about eighty. So a simple mixing 
rule cannot account for this. Anderson provides a clue to the 
explanation when he states that changing the electrode separa- 
tion by a factor of two had little effect on his readings, 
which suggests that most of Anderson’s observations were due 
to something going on at the surface of the electrode. 

Actually, this should not have been surprising. Anderson's 

collaborator, Edlefsen, reviewed an article by G.H. Cashen. 

The second point in Cashen 's summary is; 

The results (of capacitance measurements) de- 
pend on the electrodes used. With mercury, all 
soils give curves of the same general type for 
the variation of the capacity with moisture, be- 
cause the capacity effects associated with the 
soil-electrode interface are large compared to 
those due to the soil. . . The results ob- 
tained with carbon electrodes, though depending 
on the texture of the soil, generally confirm 
the. . .(results). . . with mercury electrodes. 

Apparently, Edlefsen did not fully appreciate the significance 

of Cashen 's conclusions. 
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Soon after Anderson's article appeared, however, £.C. 

Childs (9) pointed out some of these problems and stressed 
the effects of leed^y dielectrics in capacitors. Childs as- 
sumed some reasonable values of parcuneters of the model and 
showed that Anderson's results could be explained in this way. 

In 1945, E.F. Wallihan (10) of Cornell University re- 
ported on his work. He also used plaster of Paris between 
insulated electrodes. To begin, he noted that field tests 
of Fletcher's soil moisture probe were disappointing. Wal- 
lihan offered three factors which he felt were important and 
not given sufficient attention in earlier work; 

(1) The use of medium having fixed porosity as the 
material between the electrodes, 

(2) The relative merits of insulated and uninsulated 
electrodes, 

(3) Choice of frequency. 

Wallihan' reported only preliminary results cuid no subsequent 
paper of his has been found. 

In 1947, Thorne and Russell (11) published their con- 
clusions about dielectric methods of determining soil mois- 
ture. They attempted to test Child's theory directly and 
found that it was confirmed, at least qualitatively. Thus the 
apparent capacitance is a function of sample resistance among 
other things. After publication of this article cuid subsequent 
mention Of it in Advances in Agronony, few articles on di- 
electric methods appeared in soil science journals, it is worth 
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noting that Thorne and Russell used electrodes which were 
insulated with thick pieces of glass. 

Specific Aspects of Soil-Moisture Interaction 

Soil Type 

Many factors are involved in soil classification, but 
the most important is texture, which is a term referring to 
the size of particles. Table 1-1 gives a quantitative 
meaning to common textural terms. 

Real soils, of course, are seldom so homogeneous as to 
fit into only one of these categories, but most soils are 
mixtures of these types. The "textural triangle" shows how 
combinations of soil types are named (see Figure 1-1) . To 
determine the textural types, first project the silt content 
parallel to the clay side, then project the clay content 
parallel to the sand side. The region where the lines inter- 
sect is the texture type. For example, a soil composed of 
25% clay, 35% silt, and 40% sand would be classed as a loeun. 
The dashed lines on Figure 1-1 show how to project these 
various component values. The contour lines for constant silt 
content are parallel to the clay axis, the contour lines for 
constant clay content are parallel to the sand axis, and the 
contour lines for constant sand content are parallel to the 
silt axis. 

How does particle size affect the permittivity of a soil- 
moisture mixture? Water is bound onto surfaces by adsorption. 


•r 
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Table l-l. Definition of soil texture types. 


U.SD.A. System 


particle 

size 

range 

2-h 


texture type 


very coarse sand 


1-5 

5-. 25 
.25-. I 
. I-.05 
.05-.002 
<,002 
millimeters 


coarse sand 
medium sand 


fine sand 


very fine sand 
silt 


clay 


International System 


texture type 


coarse sand 


fine sand 


silt 


clay 


particle 

size 

range 


2-.2 


2-.02 


02-.002 


<.002 


millimeters 



























1 ? 



100 90 80 70 60 50 40 30 20 10 0 

% Sand 

Figure 1-1. Textural Triangle. 
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and water that is "bound" has a lower dielectric constant 
than free water in the range of frequency of interest to us. 

An example of this is the permittivity of ice, where each 
molecule is bound to its neighbors. This makes it difficult - 
for the dipolar water molecules to reorient themselves in 
response to an electric field. Since the permittivity of 
water depends largely on the motion of dipoles, the permit- 
tivity becomes very small at all frequencies except the 
lowest. The behavior of adsorbed water on clays is somewhat 
similar, according to Hoekstra (private communication). The 
permittivity of free water; ice, and adsorbed water on clays 
are illustrated in Figure 1-2. 

The amount of water adsorbed depends on the available 

surface area. As particle size decreases, the total surface 

area of a given volume of soil increases dramatically. It is 

possible to calculate the surface area per cubic centimeter 

for various size particles. This shows that surface area is 

approximately proportional to 1/radius, with a 0.0001 cm size 

2 

particle giving 31,550 cm per cubic centimeter. Actually, 
this greatly underestimates the area available since the in- 
sides of the particles also open to admit water molecules (see 
Figure 1-3). So clAys have far more surface area than sands 
and therfore adsorb much more water. Hence, a larger percen- 
tage of the water contributes very little to the dielectric 
constant, and one would expect clays to have a smaller di- 
electric constant at a given moisture content than sands or 


silts. 
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Figure 1-2. Permittivity of Water As A Function of Frequency 


a Water Clay 



Figure 1-3. Adsorbed Water 
OH; Clay. 


• Cation 

9 9 $90® 



Figure 1-4. Adsorbed Cations 
on clay. 



Figure 1-5. Polarization of Clay Particles. 
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Another problem arises from the fact that clays ad- 
sorb ions as well as water onto their surfaces (see Fig- 
ure 1-4) . This is because the clay particles have a struc- 
tural negative charge due to substitution of lower valence 
ions (magnesium, etc.) for aluminum in the crystal lattice 
and also due to broken edges in the lattice. The cations 
can move freely along the surface of the clay particle in 
response to fields. When a field is applied, therefore, 
the clay particle becomes a large dipole (Figure 1-5) . 

Hence, clays have a high dielectric constant at low fre- 

4 

quencies (generally less than 10 Hz) . Several investi- 
gators have formalized this intuitive explanation. O'Konski 
(12) introduced the notion of surface admittance to quantify 
the ability of ions to move in response to an applied elec- 
tric field. Surface admittance is merely the ratio of the 
current on the surface of a clay particle to the voltage 
across the clay particle. Schwan (13), however, showed 
that this concept was not adequate to explain experimental 
results. Schwartz (14) extended O'Konski's concept to in- 
clude complex surface admittances which imply some energy 
storage. This improved model is supported by Schwan 's ex- 
perimental results. By applying Schwartz's model to clay 
soils which have a lower limit (for particle size) of 0.2 
micrometer, it can be shown that this phenomenon has a def- 

A 

inite effect at frequencies below about 10 Hz. 
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Soil Conciition and Time 

M.L. Wiebe (15) of the Texas A & M Remote Sensing Cen- 
ter has recently completed a program of measurements on 
various soils in Texas at a frequency of 10.625 GHz. His 
results seem to indicate that variables such as homo- 
geneity, compaction, and time have a significcuit influence 
on the permittivity of soils. Table 1-2 illustrates the 
significance of various effects. A few observations about 
the factors charted in the table may serve to clarify the 
information derived from Wiebe 's graphs. Different com- 
pactions were obtained using a controlled full-compacting 
hammer. Curing time is the amount of time that the soil was 
allowed to stand in a sealed container. Only Tarrent 
Stoney Clay is not homogeneous, and the table shows the 
effect of sieving, bach entry in the table shows the 
relative influence of the different variables; the entries 
can be considered as partial derivatives evaluated at the 
midrange of moisture content. The table does show that 
these variables do have a significant influence on the per- 
mittivity. 

To get a better idea of the implications of Table 1-2, 
consider Table 1-3. This table is designed to show the 
possible inaccuracies that could result from attempting to 
infer moisture content from the permittivity of soil. Sup- 
pose one measures the dielectric constant of soil whose, 
moisture content, compaction, and curing time are unknown. 
What can be inferred about the moisture content? Table 1-3 


Table 1-2. Relative Influence of Moisture , Compaction, Curing Time 
emd Homogeneity on Soil Permittivity. 
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tne ocner symooxs. 'me ditxerence, or ch 2 mge, in compaction 1 
20N/cm^ - Siycm* and the difference in curing time is one day. 
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Table 1-3. The Possible Range of Moisture Content for litilch 

e is Constant, Provided Compaction and Curing Tine 
Can Vary Arbitrarily. 


— 
soil type 

range of moisture 

Torrent Stoney Cloy 

10% 

Gilo Sandy Loom 

2% 

Hobon Sandy Loam 

5% 

Amarillo Fine Sandy Loam 

1% 

Abilene Clay Locmn 

5% 

Houston Black Clay 

10% 

Miller Clay 

10% 

Lake Sand 

2% 


Note that this teible can be interpreted as a worst case 
error table. If more information were available, the 
error could be reduced. This additional information 
would involve the compaction and the curing time (at 
least, for laboratory samples) . Presumably, field sam- 
ples are thoroughly cured. 
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shows that the moisture content could be any number with- 
in a range of moisture contents, depending on what the 
compaction and curing time actually are. However, com- 
paction and curing time are restricted to certain limits. 

2 2 

If 20 N/cm and 5 N/cm represent extremes rn compaction 
in normal soils and variations in curing time are no 
greater than those considered here (24 hours) , it appears 
that these variables will cause considerable difficulty 
only with clay soils, as indicated by Table 1-3. 

Ion Concentration 

In a series of experiments, Mandel and Jenard (16) 
showed that salt solutions of 0.02 molarity or less dif- 
fer negligibly from water in permittivity. Beyond 0.02 M, 
the permittivity of the solution is about (80 - 10 M) , 
where M is the molarity of the solution (17) . What is the 
range of molcurities in soil solutions? Sea water is ap- 
proximately a 0.16 M solution. Clearly, the soil solution 
must have a lower molarity than this in arable regions. A 
0.16 H solution would have a dielectric constant only 2% 
less than that of distilled water. Therefore, it seems 
reasonable to expect that the presence of salt in soil 
should have little effect on the permittivity. However, 
one must be cautious with this kind of analysis since soil 
physics is very complicated and other factors may be at work. 
For this reason, it is important to use experimental data to 


25 


determine the role of various influences on the permit- 
tivity of soil. 

Temperature 


As with ion concentration, the best available exper- 
imental results deal with water, not soil and water. 
Dorsey (18) gives the following equation for the depen- 
dence of permittivity on the temperature of water, 


1-3 E = 81.47 


fT - 17| /T - 17 

1 - 4.696! /f 10.2 


1000 


1000 


1 


Where T is in degrees Celsius. 


Frequency 


The section on soil type included some information 
about the frequency dependence of the dielectric constant 

4 

of clays. Recall that below 10 Hertz, the permittivity 
of clays can be very large due to the large surface admit- 
tance of the clay particles. 

As the frequency increases above two GHz, the permit- 
tivity of water declines. According to Eisenberg (19), 
Figure 1-6 gives the frequency dependence of the dielectric 
constant of water. 
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Conclusions 


The conclusions from both of the previous sections 
are summarized below under their subheadings. 

Determination of Soil Moisture Using Dielectric Methods 

From the experiments of Anderson and Thorne and Russell 
it seems very probable that electrode polarization is re- 
sponsible for many of the observed effects. This implies 
that resistance veuriations of the soil with water have a 
strong influence on what is actually measured. This is a 
serious limitation since resistance is a function of ion con- 
centration as well as of water. The concentration of ions 
can vary considerably, even in one location. Fertilizer 
treatments, for example, can have a large effect on the 
measured resistance of the soil. 

Fletcher claimed that salt had little effect on his re- 
sults. However, Fletcher used a plaster of Paris block be- 
tween the electrodes. The plaster of Paris is soluble. In 
fact, the block is so soluble that it dominates the situation 
and determines the ion concentration between the electrodes. 
Therefore, variations in salt content have little effect. 
Since it is the resistance, or something proportional to it, 
that Fletcher’s method measures, his method does not seem to 
have any real advantage over resistance techniques. 

In conclusion, most of the investigators cited either 
failed to take electrode polarization into account or to elim 
inate its effects. 
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Specific Aspects of Soil Moisture Interaction 

The work of Mandel and Jenard, and others, seem to 
indicate that ion concentration may not be an important 
factor once electrode polarization is eliminated. 

According to the literature, soil type should have 
an influence on the permittivity. For a given moisture 
content, clays should have a lower permittivity. 

The best frequency range in which to make measure- 

4 9 

ments is probably between 10 Hz and 10 Hz. Below this 
range clays can have a very high permittivity which is 
not related to their moisture content. Above this range 
the dielectric constant of water begins to fall off. As 
the permittivity of water decreases, the permittivity of 
the mixture is influenced more by the properties of the 
soil. At a few megahertz water has a dielectric constant 
20 to 40 times larger than that of most soils whereas at 
20 GHz the permittivity of water is only 5 to 10 times as 
large as that of most soils. Hence the “noise" increases 
considerably in this region. 
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CHAPTER 3 

PRELIMINARY INVESTIGATION OF 
THE PERMITTIVITY OF SOILS 


Early in the course of this project , it was decided 
to measure the permittivity of some soils to obtain data 
not available in the literature. The program of measure- 
ments explained a great deal about electrode polarization 
and made possible the interpretation of many of the pre- 
vious investigations. 

Apparatus and Measurement Procedure 

A Wayne-Kerr B601 bridge together with its associated 
detectors and oscillators was used in this investigation. 

A circuit diagram is shown in Figure 1-7 and a block dia- 
gram in Figure 1-8. 

The operation of the bridge was quite simple , consisting 
of an initial balance with the sample disconnected and a 
final balance with the sample in the circuit. The capaci- 
tance and resistance of the unknown (containing the soil as 
a dielectric) could be read directly from the balance knobs. 

The soil was placed in a scunple holder whose capacitance 
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Figure 1-7. Circuit Diagram (0.02 - 2 MHz range) 
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and resistance were measured by the bridge. The sample 
holder underwent modification during the course of the in- 
vestigation. In the beginning brass was used but it tarn- 
ished rapidly. Then a nickel-plated holder was tried which 
also reacted with the sample. The sample holder illustrated 
in Figure 1-9 was constructed of stainless steel. It seemed 
to resist chemical reaction and gave results which agreed 
with those obtained by using platinum electrodes. The 
sample holder is a coaxial capacitor with guard rings which, 
serve to eliminate fringing of the electric field of the ca- 
pacitor. In the absence of electrode polarization, the per- 
mittivity of the sample can be calculated by dividing the 
measured capacitance by the capacitance of the sample holder 
when the soil is removed leaving air as the dielectric. 

Preparation of Samples 

In the beginning, preparing soil scunples presented a 
real problem. Samples with uniform moisture content were 
needed and, in addition, we wanted to be able to mix exactly 
the desired moisture content. Mixing water and soil directly 
does not give a uniform mixture. After much trial-and-error 
the following procedure was chosen. Soil was dried in an 
oven and placed in a freezer after being weighed. Per- 
centage moisture content is defined as (weight of water)/ 
(weight of dry soil) x 100. Using this definition, the amount 
of water required was calculated then increased by 5% to cover 
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Figure 1-9. Soil Sample Holder (0.02 - 2 MHz range} 
Center conductor is 19.05 cm long. 
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losses. This percentage was determined by experience. The 
water was measured with a syringe and placed in five-cc 
aluminvun foil containers. The water was taken to the free- 
zer and frozen. After several days, the ice was crushed 
and scraped into a jar where it was thoroughly mixed with 
the dry chilled soil. Since the freezer was maintained at 
-15* to -20"P, very little melting took place. The mixture 
was then poured into the sample holder. This was removed 
from the freezer and allowed to come to equilibrium for 12 
to 24 hours. 

Results and Discussion 

The data obtained in this investigation are tabulated 
in Appendix 1. Some representative data are shown in Fig- 
ures 1-10» 1-11, and 1-12. Several; points are worth noting. 

f ' ■ ‘ 

First, the capacitance is definitely a function of frequency. 
Second, the values one might infer for the permittivity of 
the soil-water mixture turn out to be considerably more than 
eighty in a number of cases. In addition, the addition of 
salt, even in small amounts, resulted in an increase in the 
capacitance of the sample. These observations led to the con- 
clusion that electrode polarization was a significant factor 
in the results. This point will be covered in more detail 
in the next chapter. 

In this series of measurements the values of capacitance 
occasionally varied with time. Repeated measurements on 

1 
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Figure 1-12. Capacitance Vs. Frequency, Sand, 5.8% Moisture Content 
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the empty scunple holder and on commercial capacitors and 
resistors proved to be completely repeatable and drift free. 
Eventually it was concluded that the time- varying values ' 
were partly due to chemical reactions at the electrodes 
and partly due to something happening within the sample. 

The observations are summarized below. 

Observations 

(1) On muck soils in a brass holder the capacitance 
increased with time, sometimes by as much as 10% in the 
course of an hour. Some corrosion of the electrode was 
evident, but it was not possible to correlate the amount 
of corrosion with the time duration. 

(2) A similar s^unpje holder was nickel plated. The 
capacitance of Bentonite sample 5 (Bentonite is a clay min- 
eral) decreased about 10% in one hour. A green residue 
was noted on the electrodes. (It may have been NKOH)^ or 
some similar nickel compound. Many of them are green). 

In order to look into this a little more deeply, two ident- 
ical samples of Bentonite wetted to the wilting point were 
prepared and the following experiments were performed. 

Sample 1, Experiment 1 

Connect the sample holder to the bridge and balance 
the bridge. Changing nothing else, simply re- 
balance the bridge every five minutes. 
sample 2, Experiment 2 

Connect t)ie sample holder to the bridge and balance 
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the bridge. Rebalance the bridge every five 
minutes but disconnect the power source from 
the bridge between measurements. 

Sample 1, Experiment 3 

Seunple 1 was allowed to sit for two hours fol- 
lowing Experiment 1. Attach it to the bridge 
and rebalance as often as possible. 

Sample 2 , Experiment 4 

Attach Sample 2 to the bridge and balance the 
bridge. Rebalance the bridge every five minutes / 
disconnecting the sample from the bridge between 
measurements. 

The data obtained from the first three experiments are shown 
in Figure 1-13. In Experiment 4, the capacitance varied 
little from measurement to measurement. A voltmeter showed 
a very small voltage across the terminals of the sample holder. 

(3) A stainless steel sample holder was made. The cap- 
acitance of a Bentonite sample decreased 3% in an hour. No 
corrosion at all was observed on the electrodes. 

(4} The capacitance of silts and sands without much or- 
ganic matter did not vary with time. 

(5) The capacitance of Bentonite samples measured using 
platinum electrodes decreased about 3% in an hour. 

Recall that in each case where the sample holder had an 
electrical discharge path (through the bridge) , the capacitance 
values changed with the passage of time. This fact, together 
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Figure 1-13 Change in Capacitance With Time, of Soil-Filled Test Container 
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with the forming of a residue on the electrodes and the 
existence of a voltage between the electrodes led to the 
conclusion that electrochemical reactions at the electrodes 
were responsible for the observed time variation with brass 
and nickel sample holders. Platinum, however, is reputedly 
quite inert. Both platinum and stainless steel sample holders 
gave the same time varying results when they were filled with 
clays. These results were independent of connection with the 
bridge. Possible explanations for this phenomenon will be 
discussed in Chapter 6. 

The next chapter deals with electrode polarization. An 
understanding of this concept aided in interpreting the re- 
sults of the preliminary investigation and the planning of sub 
sequent investigations. 
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CHAPTER 4 

ELECTRODE POLARIZATION 


Introduction 

Due to the work of many electrocheinists, electrode 
polarization is now fairly well understood. Many details 
are not clear, however, and some concepts are controversial. 
Basically, the explanation is as follows. 

When an electrode is inserted into a solution a potential 
is set up between the electrode and some of the ions in the 
solution. These ions are attracted to the electrode, but, 
owing to their finite radius, can approach no closer than a 
few angstroms. These ions, together with the charge they in- 
duce in the electrodes, constitute the "diffuse double layer." 
Maxwell's laws can be applied to this system, taking into 
account the facts that the charge is discrete, the ion layer 
is diffuse, the dielectric constant of the solution varies, etc. 
When all the parameters are known or can be reliably estimated, 
the analysis yields results that are in fair agreement with 
experiment. This kind of agreement has been achieved only for 
pure solutions where impurities are constantly removed via 
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a mercury'drop electrode. The addition of just a little 
impurity (for example, a small cunount of gelatin (20) } 
changes things considerably. One of the best references 
for the detailed theory of electrode polarization is Bock- 
ris and Reddy (17). 

A Model for Electrode Polarization 

In view of the fact that this model is so difficult to 
handle cUialytically and of dubious value for the kind of so- 
lutions used in this work, a somewhat simpler model was de- 
veloped, as discussed in the following, ions are the main 
charge carriers in soil. These cannot penetrate the metal 
electrode so they accumulate at the surface of the electrode. 
There is also some conductance. The model for these effects 
is a capacitor in parallel with a resistor. So the sample 
holder containing soil is modelled by the circuit shown in 
Figure 1-14. The resistors and capacitors in this model are 
not fixed*fvalue lumped elements. For the moment, let us as- 
sume that and are constants. Intuitively, we expect 
and R^ to depend on frequency since these parameters de- 
pend on the motion of ions, which have considerable inertia. 
Schwan (21) used the model shown in the lower half of Figure 
1-15. He claimed that the frequency dependence is as shown 
in the second pair of equations in Figure 1-15. The first 
pair of equations relates the constants of the two circuits 
of Figure 1-15. We can check Schwan' s estimate of by 


44 


1— A/VV — 


fAAAn 

Ri 


Ra 




— 


C| Cg 


1 R,+ l/jo>C| s electrode sheath impedance 

2 R^+ l/jojC^- sample impedance 
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making an estl^te of our own. of Figures 1-14 and 1-15 

can be thought of as the capacitance of a parallel-plate ca- 
pacitor. One plate is the metal electrode while the other is 
the ion sheath which covers the electrode. So » eA/d, A 

being the area of each plate. We will calculate the capaci- 

2 

tance on a unit area basis, thus A = 1 cm . The separation 
distance is d. This is about twice the hydrated radius of 
the ions, or about 8 angstroms. The dielectric is water, but 
the water is bound, and this reduces its relative dielectric 
constant to around 10 (17) . So = 11 yF. (See Figure 1-16 
for this detailed calculation.) A detailed theoretical cal- 
culation (17) yields » 16yF, independent of ion type and, 
within certain limits, of ion concentration. Now that the 
parameters are better understood, the equivalent circuit for 
electrode polarization may be analyzed. The detailed analysis 
is in Appendix ill. The actual values measured are R and 
the eqpiivalent parallel components of the circuit of Figure 
1-14. Expressions for C and R are given in Figure 1-17. In 
Figure 1-18 R and C are sketched as functions of frequency. 

This relationship can be tested by substituting into the 
equations for C and R reasonable values of R^ and C 2 and plot- 
ting values of C and R versus u. These should agree with the 
experimental values which have been given in Chapter 3, 

Looking at Figure 1-10, R^ is about 260 n and about 60 pF. 

2 

The aurea of the small electrode is about 10 cm , so is 

-4-7 

approximately 2 x 10 f * F. Substituting these values into 
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Figure 1-17. Equivalent Parallel R and C For Circuit of 

Figure 1-14 (Sample Holder Containing Soil). 



Figure 1-18. R and C of Figure 1-17 as Functions of Frequency. 






the equation for C in Figure 1-17, the values plotted on 
Figure 1-19 are obtained. The agreement is fairly good, con- 
sidering that is calculated from basic considerations. 

Were a little larger and slightly less dependent on fre- 
quency, the agreement would be nearly perfect. 

Conclusions 

\ 

The formula for C derived here allows one to estimate the 
effect of ion concentration in the soil sample. This concen- 
tration strongly influences the conductivity of a seunple but 
should have little effect on its permittivity. Thus, the sam- 
ple resistance is a function of ion concentration. As the 
equations of Figure 1-17 show, the overall capacitance, C, of 
the soil- filled sample holder is a function of the resistance 
of the soil inside. It appears from the equation for C in Fig- 
ure 1-17, that each decade drop in the resistance, R 2 » shifts 
the C vs. f curve about 1 1/2 decades in frequency. 
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CHAPTER 5 

INVESTIGATION OF THE PROPERTIES OF 
SOME SOILS 


In order to eliminate the effects of electrode polar- 
ization, measurements were obtained at frequencies greater 
than 2 MHz. In the beginning, a range of 5 to 40 MHz was 
used, but later it was necessary to use still higher fre- 
quencies. 

Measurements in the 5-40 MHz Range 
Apparatus and Measurement Procedure 

The instrument used to perform the measurements was a 
1606-B General Radio R-F Bridge. In addition to the bridge 
itself, various oscillators and detectors were required. A 
schematic diagram is shown in Figure 1-20. 

The soil was placed in a coaxial capacitor, shown in Fig- 
ure 1-21. The capacitor was made of nonmagnetic stainless 
steel, machined to the same dimensions as a section of General 
Radio 50-C transmission line: outer conductor, 1.587 cm OD, 

1.429 cm ID; inner conductor, 0.621 cm DIA. The coaxial ca- 
pacitor was 12.7 cm in length. At one end of the capacitor a 
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Figure 1-20. Block Diagram (5 - 40 MHz Range) 
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Figure 1-21. Sample Holder (5 - 40MHz Range) 
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General Radio 874 Locking Connector was attached, enabling 
the sample holder to be attached directly to the GR 1606-B 
R-F Bridge. 

The procedure for performing the measurements was adapted 
from the manual for the 1606-B. The instrument has four ad- 
justable dials; two are initial balance controls, and the 
other two indicate sample resistance and reactance. (The re- 
lationship between the measured reactance and the permittivity 
is stated in the next section.) Measurements were taken at 
5 MHz, 10 MHz, 20 MHz, and 40 MHz. After the last measurement, 
the sample holder was disassembled and the contents emptied in- 
to a small container. By using the gravimetric method (see (22) 
for details of this method) the moisture content of the seimple 
was obtained. 


Calculations 

The instrument measures the reactance and the resistance 
of the sample holder. The permittivity and conductivity can 
be calculated by using these quantities, and in this section the 
equations used to calculate the permittivity are set forth. Also, 
a test case using water was worked to show that the equations do 
indeed yield the permittivity. 

Admittance. The first step is to convert the series re- 
sistance, R, and reactance, X, (impedance) to the equivalent 
parallel capacitance, C, and conductance, G, (admittance). The 
measured quantities are R and Xf where f is the frequency in 
MHz. The desired quantities are C and G. Let D denote -Xf. 
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1-4 


lr5 


R - -1 — = . 

J f G + jwC 

^ -D i Q ^ Rf^ 

+ D* 2tt 10®/ 


In the above equations/ the unit of R and X is ohms, that 
of C is farads, and that of G is mhos. 

Fringing Capacitance Correction. In this Scimple holder 
there is a small amount of fringing capacitance which must be 
subtracted. A computer program was written to calculate the 
fringing capacitance, which was found to be independent of the 
dielectric constant and equal to about 0.38 pF. This value is 
simply subtracted from the calculated capacitance. Further 
details on this calculation are given in the second part of 
this chapter. 

Permittivity Calculation. Next, it is necessary to find 
a relationship between permittivity and capacitance. If the 
capacitor could be treated simply as a lumped capacitor, per- 
mittivity could be calculated simply by dividing the measured 
capacitance by the capacitance of the empty sample holder. 

Even though this gives good results, for a more exact calcula- 
tion it is necessary to ta}ce into account the distributed 
nature of the capacitance. The sample holder can be treated as 
a section of transmission line terminated in an open circuit. 
The input admittcuice of the sample holder is 


Y_ + Yf, tanh yl 

Y = Y ^ ^ 

IN 0 Yq + tanh yl 


where 


= input admittance, 

Yn = characteristic admittance. 


DwC. + G( 


= load admittance) for open circuit Yj^ = 0, 

= propagation constant, /(jwCj^ + Gj^) (jwLj^ + f 
capacitance per unit length, 

= inductance per unit length, 

= resistance per unit length, 

= conductance per unit length . 

For an open circuit load, 

1-7 

If tanh yH = yl ,then 


Y 

R, 


^IN = 


1-8 


Yin “ 


One might reasonably ask whether yZ is sufficiently small 
to warrant the assumption that tanh y^I - Y^^* Recall that 
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Y = VfjuL + R ) (jwC. + G.) , 


Z £' — I 

is determined by geometry and permittivity; 

1-10 C. = . 

Z Zn b/a 


where 


e = permittivity, 
b = dia. of outer conductor, 
a = dia. of inner conductor. 

Evaluating C, 

1-11 C = ZCj^ = 0.667e^pF. 

We know that for an empty sample holder, 

1-12 Y = j“ /lC = j (j/3 X 10®, 

(e^ = 1 for an empty sample holder) 


I 
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1-13 LC = , 

9 X 10^^ 

1-14 L = 1.66 X 10 ^ Henry. 

is determined by the conductivity of the soil and geometry. 

' G = 7.57<,. 

R is the series resistance of the transmission line. Since the 
sample holder is made of stainless steel there may be some rea- 
son to doubt that R is small. The resistance per unit length 
(R^) can be calculated by using the skin depth approximation. 
See Figure 1-22. 



where 


p = the resistivity of the material, 
A = the cross sectional area. 

1-17 A = 62ira 

where 


a 


1-18 

where 


the diameter of the inner conductor, 


6 




0.063 meter. 


0 ) = the frequency in radians/second, 

“7 

U = the magnetic permeability, 4ttx10 Henry meter, 
p = resistivity, 0.72 x 10 ^ ohm-meter. 

So Rj^ = 4.25 X lO"^ /f, 
where 


f = frequency in Hertz 
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a = radius of solid conductor 
S = skin depth 

A = cross-sectional orea supporting conduction 

Figure 1-22. Skin Depth of Conductor. 
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Comparing to 

1-19 = f X 1.045 X 10"S 

1-20 = /T X 4.25 X lO"®. 

Hence, = 1.0 while = 0.0042 at approximately 1 MHz. 

It can be concluded, therefore, that R is negligible 
compared to coL. Thus can be expressed as 

1-21 yl - (jtuC + G). 

Note that this is a complex number. Hence, when calculating 
tanh yi, the formula for complex arguments must be used. 


1 n-, _ sinhta) + jsin(b) 

1-22 tanh (a + jb) - cosh (a) + jcos(b) - 

Using this formula, tanh be expressed as a series ex- 

pansion in terms of the real and imaginary parts ot yl. 
Finally, can be expressed as a series expansion. Keeping 
two terms , 



Measurements on Water. As a test of the correction formulae, 
measurements were performed on the sample holder using water as 
the dielectric. The results are summarized in Table 1-4. 
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Table 1-4. Measured Permittivity of Water. 
Distilled Water, Conductivity = 0.003 Mho/meter 
Frequency in MHz Relative Permittivity 

5 80.2 

10 79.5 

20 78.0 

40 75.5 

0.005% NaCl Solution, Conductivity = 0.02 Mho/meter 
5 79.5 

10 77.0 

20 78.0 

40 76.5 

0.01% NaCl Solution, Conductivity = 0.03 Mho/meter 
5 81.3 

10 81.3 

20 81.5 

40 75.0 

0.02% NaCl Solution, Conductivity = 0.05 Mho/meter 
5 80.8 

10 82.3 

20 83.0 

40 75.0 

Q.1% NaCl Solution^ Conductivity =0.2 Mho/raeter 
5 88.5 

10 80.4 

20 93.0 


40 


88.2 
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Results 

The complete results for the capacitance, resistance, 
permittivity, and conductivity of several different soils, 
for various moisture contents, and for the entire range of 
frequencies, are tcibulated in Appendix I. Figures 1-23, 1-24, 
and 1-25 illustrate some representative data. The samples 
used to obtain the data were prepared according to the pro- 
cedure outlined in Chapter 3, Due to the smaller size of 
the sample, only four to eight hours were required for the 
samples to come to thermal equilibrium. These results are 
discussed in Chapter 6. 

Measurements in the 250 - 450 MHz Range 

After the measurements on the laboratory-prepared samples 
were completed, it seemed logical to attempt to measure un- 
disturbed samples from the field. Unfortunately, it was not 
practical to do this using the equipment for the 5-40 MHz 
range. In order to obtain the required minimum capacitance, it 
was necessary to use a sample holder which was 12.7 cm in length 
Soil samples of this length could not be forced into the sample 
holder without binding, compaction, and unreliable contact. 

Good results could be obtained, however, for soil samples 
of one centimeter in length. Samples of this length were suit- 
able for use in the 250 - 450 MHz range. In addition, it would 
be helpful to have information on the permittivity in this 
higher frequency range. For these reasons, it was decided to 


39.8 



Figure 1-23. Log e^. Vs. Frequency for Three Miami Silt Loam Samples 



Figure 1-24. Log Vs. Moisture Content, 0.2% Salt Solution (Miami Silt Loam) 




Figure 1-25. Effect of Salt on e of Miami Silt Loam 
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perform measurements on laboratory-prepared samples and field 
samples in the frequency range of 250 -450 MHz. 

, Apparatus and Measurement Procedure 

A General Radio 1602-B Admittance Meter was used in this 
investigation. A block diagram of the instrument, together 
with its associated oscillators ^lnd detectors is shown in Fig- 
ure 1-26. The Admittance Meter has a four-port junction at its 
center. Connected to the appropriate ports are an oscillator, 
a known conductance, a known susceptance, and the unknown. Three 
loops inside the junction are connected to the detector, which 
can be nulled by varying the coupling of these loops. When a 
null has been achieved, the coupling of the loops is related to 

the input admittance of the unknown port in a known fashion. 

The Admittance Meter is a particularly easy instrument to 
use. The following procedure was adapted from the manual for 
the instrument. To begin, the source oscillator is set to the 
desired value, and the auxiliary oscillator is set to obtain 
maximum deflection on the detector. Occasionally, it is nec- 
essary to observe the mixer current. If the mixer current is 
too low the coupling loop on the auxiliary oscillator can be 
adjusted to give better results. After the susceptance stub 

is set according to the frequency, then the seunple holder is at- 

tached to the "unknown" port. The three levers attached to the 
coupling loops are adjusted to obtain a null in the detector, 
and the positions of the levers are recorded. 

Several points deseirve attention. First, it is important 
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Figure 1-26. Block Diagram (250 - 450 MHz Range) 



.58cm 00 tubing Drilled and tapped 
(303 stainless steelh to occept GR 874 connector, 

hic"^ including grooves ond slot shown 




0.62 cm 01 A Rod 


Figure 1-27. Sample Holder (250 - 450 MHz Range). 
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that the source oscillator be set as accurately as possible, and 
that the frequency of the oscillator be checked from time to 
time by means of a slotted line, unless the oscillator is known 
to be highly stable in frequency. Small errors in frequency 
cause errors in the final result, Secondly, the susceptance 
stub should be set according to the marks on the stub. The 
procedure given in the GR manual for setting the stub when the 
frequency is unknown leads to errors and should not be used. 


Calculations 


Five sets of calculations are performed on the data to ob- 
tain the permittivity: a line length correction, a finite 

sample length correction, a fringing capacitance correction, an 
optional electrode polarization correction, and a calculation 
of the permittivity from the seimple input admitt 2 mce. 

Line Length Correction. The Admittance Meter measures the 
input admittance at the junction, but the sample is approximately 
6 cm away from the junction. The equation which relates the 
input admittance at one point in a lossless line to the input 


admittance at another point is 


1-24 


V = V 

'‘in 1 *0 - ] Y, 


where 


Yjj^ ^ = input admittance at the desired point, 
Yjjj 2 “ input admittance at the known point, 

Yq = characteristic admittance of the line, 
3 = phase constant, . 

= distance between point 1 and point 2. 


£ 


65 


Note that Yq and 6 are characteristics of the transmission line 
and, therefore, known. 2 measured by the Admittance 

Meter. 

In order to use the formula it is necessary to determine 
JL, which is the distance between the junction and the sample. 
This can be accomplished by attaching a known admittance to 
the "unknown" port. The known admittances used in this inves- 
tigation were GR W03 (open circuit) and GR WN3 (short circuit) 
terminations. The use of these terminations not only simpli- 
fies the formulae but also provides a check on the error. 

The simplified formulae for line length corrections are 

(short circuit) Y^^^ ^ - Y^/tan Bt, 

(open circuit) 2 ~ ^0 

Note that the product of the measured input admittances for 

2 

short circuit and open circuit is just Yq. Table 1-5 shows the 
data collected using WN3 and W03 terminations as well as the 
calculated, value of I and the estimated error. Knowing i, one 
can calculate the input admittance at the seunple holder by using 
Equation 1-24. 

Finite Length of the Sample Holder. In the section on mea- 
surements in the 5-40 MHz range, equations were derived to re- 
late the input admittance to the capacitance and conductemce 
per unit length of the sample holder filled with soil. The re- 
lationships were 

Im [Y^j^l 

'i + 

3 

i J 


1-25 


C 
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Re [Y^^] 

1-26 G = 


fl + ^ 

fc - ^]] 


i J 4 


These relations are also useful in the present calculations. 

Fringing Capacitance. The correction for fringing capaci- 
tance mentioned in the section on measurements in the 5-40 
MHz range is explained in more detail here. To begin, let us 
consider what fringing capacitance is and how it might affect 
the results. Looking at the s^ple holder, one can see that 
it is a coaxial capacitor. Because the capacitance per unit 
length of coaxial capacitors is well known, and since the length 
of the sample holder is also known, one might be tempted to 
conclude that the capacitance of the sample holder is just its 
length times the capacitance per unit length. This, however, 
is not the case. In fact, the capacitance of the sample holder 
is greater and the excess is termed the fringing capacitance. 

The reason for the discrepauicy is that the formula for the cap- 
acitance per unit length was derived for an infinitely long 
capacitor. Since the sample holder is finite in length the uni- 
formity of the infinite capacitor is gone. At the edge of the 
capacitor there are fringing fields which give rise to the 
fringing capacitance. See Figure 1-28. The calculation of the 
fringing capacitance consists of calculating the capacitance of 
the entire sample holder and comparing it to the length of the 
sample holder multiplied by the capacitance per unit length cal- 
culated for the case of an infinite capacitor. The calculation 
of the actual capacitance is accomplished by using the following 
relationships. 
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1-27 


C = -2 
A4> 


where 

C = capacitance, 
q = charge on the electrode, 

= voltage difference between electrodes. 
Gauss' law yields 
1-28 eE • 3s = q 


where 

— >■ 

E = electric field intensity, 

= incremental surface area, 

S = any surface which encloses q. 

Recall that 
1-29 E = -V4». 

Finally, satisfies the Laplace equation in charge-free regions 


1-30 = 0. 

Using a computer one can solve the Laplace equation by 
means of finite difference methods. (See, for example, Ramo, 
Whinnery, and Van Duzer (23).) The Laplace equation becomes a 
difference equation which can be expressed as 


1-31 

ir j 


■ 1 +(l) . . . . , 

^ ^1+1,] ^1-1,3 ^3-,]+l 


This difference equation is solved for different assumed values 
of permittivity in the sample holder. The potential found in 
this way is substituted into Equation 1-29 and the fringing 
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capacitance calculated. The result found was that the fringing 
capacitance was equal to about 0.38 pF, independent of the per- 
mittivity of the soil in the sample holder. One of the pro- 
grams used to calculate the fringing capacitance is given in 
Appendix 4. 

Electrode Polarization Correction. This correction can be 
considered optional since it produced only small changes in the 
data. Recall the formula for measured capacitance derived in 
Chapter 4 and given in Figure 1-17 » 




Calculation of the Permittivity. Once the capacitance of 
the sample holder filled with soil has been found, the permit- 
tivity can be found simply by dividing the sample holder cap- 
acitance by the capacitance in air, 0.667 pP, 

1-34 = C 2 / 0.667 x lo"^^ 

Test of Formulae. To insure the correctness of the for- 
mulae used to calculate the permittivity and to correct for 
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various effects, measurements were performed on materials 
with known dielectric constants: water, methyl alcohol, and 

plexiglas. The results are listed in Table 1-6. A computer 
program, written to obtain the permittivity can be found in 
Appendix IV. 

Table 1-6 Measured Permittivity of Water, Methanol, and 
Methyl Methacrylate. 

Relative peimittivity of water (250 MHz, 20“ C) 

Sample Holder Measured Generally Accepted 

Permittivity Value for Permittivity 


#1 80.5 80.3 

#2 79.6 80.3 

#3 80.5 80.3 


Relative permittivity of methanol (250 MHz, 20° C, Sample Holder #1) 

Measured Generally Accepted 

Permittivity Value for Permittivity 


32.2 

Relative permittivity 

32.2 

of methyl methacrylate (Sample Holder #1) 

Frequency in Maz 

Measured 

Permittivity 

Generally Accepted Range 
Values for Permittivity 

250 

2.58 

2.2 - 3.2 

300 

2.52 

2.2 - 3.2 

350 

2,41 

2.2 - 3.2 


450 


2.39 


2.2 


3.2 
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Sample Preparation 

In this investigation three different kinds of soil 
samples were prepared: laboratory-prepared samples, laboratory- 

prepared samples with salt, and seunples from the field. The 
laboratory-prepared samples were prepared according to the 
procedure outlined in Chapter 3. when salt was required it- 
was added to the soil before the ice was mixed with the soil. 

The salt content is expressed as a percent soil solution, e.g., 
a twelve gram soil sample with three grams of water added and 
0.03 gram of salt would be a sample with 25% moisture content, 
auid a 1% salt soil solution. In other words, the soil is wet- 
ted with a 1% salt solution. 

As the name implies, field samples were taken directly 
from the field, specifically, three pits dug on the Purdue Uni- 
versity Horticulture Farm. To obtain a sample, one removed ap- 
proximately three inches of soil from one face of the pit and 
inserted the outer conductor of the sample holder into , the side 
of the pit. Then the outer conductor, together with the soil 
inside it, was removed, packed in water-proof bags, eind taken 
immediately to the laboratory. In the laboratory a hole was 
drilled in the center of the plug of soil, cuid the center con- 
ductor was carefully inserted, then tamped in to insure good con 
tact with the electrodes. Finally, a special tool was inserted 
to remove excess soil so that the soil sample was precisely one 
cm in length. A rubber stopper was placed on the end of the 
outer conductor to prevent evaporation. Then, the measurements 
were performed immediately, using the Admittance Meter. 
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Unfortunately, the very wet conditions which prevailed 
during the investigation forced modification of the procedure. 

One pit remained flooded for the entire season and was not used. 
Samples from the other two pits were at or near field capacity 
of moisture. Although these data were useful, the decision 
was made to test a greater range of moisture content. The 
procedure, therefore, was modified. 

I 

A three-inch diameter brass pipe was inserted a few inches 
into the face of the pit. This was taken to the Laboratory where 
the soil plug was removed and sliced into 1.5 cm thick slabs. 
These were placed on a pressure plate and dried to the appro- 
priate moisture content. (This is a standard technique. See» 
for example, (24). ^ Once the slab was at the desired moisture 
content, the outer conductor was inserted into it, and the pro- 
cedure thereafter was identical to the one explained previously. 

The results are tabulated in the third section of Appendix 

I. Figures 1-29 to 1-34 illustrate the relationship between 
moisture content and the logarithm of permittivity. Figure 1-35 
shows the frequency dependence of the permittivity. Some pro- 
perties of the soils used in this investigation are summarized 
in Appendix II. 

RESULTS 

Several observations can be made about the data. 

1. The log of the permittivity increased approximately 
linearly with increasing moisture content. 

2. The salt content does not seem to affect the results. 

3. In this frequency range the permittivity has very 
little dependence on frequency, although the clay soil 
seems to exhibit more dependency than the others. 



74 > 



Content for Miami Silt Loam, Laboratory Prepared 
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CHAPTER 6 

DISCUSSION OF RESULTS 


In this chapter results from the investigation of permit- 
tivity of some soils in three frequency ranges are discussed. 

Measurements in the 0.02 - 2 MHz Range 

The procedures used to obtain the results discussed in 
this section can be found in Chapter 3. The data are tab- 
ulated in Appendix I. The apparent dielectric constants can 
be obtained from the capacitance by dividing the capacitance 
listed in the table by the capacitance of the empty sample 
holder, 1.65 pF. This has not been done because the results 
are not of practical value to the present investigation. 

Several aspects of the data are worth noting. The ap- 
parent dielectric constants are very high for almost all the 
soils tested. In addition, the apparent dielectric constants 
are strongly dependent on frequency. Furthermore, the conduc- 
tivity also has an effect on the apparent dielectric constant. 
These phenomena are not consistent with the assumptions that 
the apparent dielectric constant is, in fact, the actual 
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dielectric constant, and that the dielectric constant of a 
soil sample can be calculated from the dielectric constants of 
the cons^tituent soil and water by some simple mixing rule, e.g., 
e = (fraction of water)+ e (fraction of soil). 

J. SOXx 

In Chapter 4 it was shown that the theory of electrode polar- 
ization not only qualitatively explains the phenomena, but also 
makes quantitative predictions which are close to experimental 
results. 

Because the model of electrode polarization gives results 
which agree reasonably well with experimental results, one 
might ask if the theory could be used to correct for electrode 
polarization effects in the data from the first part of Appendix 
I. Unfortunately, as shown below, the answer is no. An ap- 
proximate correction formula has been derived for the case where 
electrode polarization is a minor effort. (See Chapter 5, 
Equation 1-33.) 


'AC 


'AP 




where 

e = actual dielectric constant 

= apparent dielectric constant (measured) 
w = angular frequency 

= resistance of sample in ohms 

= electrode capacitance in farads, due to an ion sheath 
on the electrode 

Cq = capacitance in farads, of empty sample holder 
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If were equal to lO'^F, equal to 1.7 x 10”^ ®F, R 2 
equal to 1280 w equal to 2it x 10® radians/second, and 
the actual dielectric constant equal to ten, the above 
formula would predict an apparent dielectric constant of 
100. Suppose the dielectric constant actually measured 
was 110. This is within 10% of the predicted value, which 
is reasonably good agreement. If one calculates the 
by using the measured value E^p = 110, one obtains 20 which 
is 100% greater than the actual value. The basic problem 
here is that the difference between two large nearly equal 
numbers always has fewer significant digits than the two 
original numbers. 

Estimating the electrode capacitance is very difficult, 
since it is proportional to the surface area of the elec- 
trode in contact with the sample. Due to surface irreg- 
ularities, the area of the electrode can be an order of mag- 
nitude larger than a perfectly smooth surface of the same 
dimensions. Even neglecting the problem of estimating the 
surface area of the electrode, one is still faced with the 
problem of estimating the amount of surface area which is 
actually in contact with the soil sairiples. Although care- 
ful packing of the soil is helpful, some variation of the 
actual contact area is unavoidable. Finally, the electrode 
capacitance itself has some frequency dependence which com- 
plicates the situation. The conclusion is that, where elec- 
trode polarization has a large influence on the data ob- 
tained, it is probably not feasible to correct for electrode 
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polarization. At higher frequencies, where and are 
not greatly different, correction is feasible but not par- 
ticularly helpful. 

An interesting phenomenon encountered in the course of 
this investigation was the time dependence of the data, de- 
scribed in Chapter 3. Even after the possibility of a chem- 
ical reaction at the electrode was eliminated, the capacitance 

4 

of clay and muck soil decreased with time. One possible ex- 
planation is that some of the moisture is adsorbed into the 
surface of the soil particles and, therefore, the permittivity 
of the sample is reduced since adsorbed water has a lower di- 
electric constant. Another possibility is that as moisture 
is adsorbed into the surfaces, less moisture is available to 
support ionic conduction. This would result in a decrease in 
sample conductivity with a consequent decrease in the effect 
of electrode polarization. Hence, the apparent dielectric con- 
stant would decline. 

Actually, these two explanations are not mutually ex- 
clusive and, in fact, there is some evidence that both pro- 
cesses occur. Wiebe (15) , whose investigation was cited in 
Chapter 2, noted a decline in measured dielectric constant with 
time. His measurements were performed at a frequency high 
enough to preclude any electrode polarization effects. On the 
other hand, decrease in conductivity of laboratory seimples 
with time has been noted by other authors (25) , and the model 
in Chapter 4 shows that this would decrease the apparent per- 
mittivity. This effect should not be present, of course, with 
soil in the field. 
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This preliminary study served to explain several studies 
published in the literature. For example, the very high ap- 
parent dielectric constants measured by Anderson (8) must 
have been the result of electrode polarization. Likewise, the 
strong influence of salt content noted by Thorne and Russell 
(11) Ccin be traced to electrode polarization, which illustrates 
the futility of insulating electrodes, as Wallihan and Thorne 
and Russell (10, 11) did in order to overcome the effects of 
soil conductivity. The insulation merely reduces the elec- 
trode capacitance somewhat, which increases the range of fre- 
quency over which electrode polarization has an important ef- 
fect. 

Measurements in the 5 - 40 MHz Range 

The procedures used to obtain the data discussed in 
this section are given in Chapter 5. The data are tabulated 
in the second section of Appendix I. 

In an attempt to overcome the effects of electrode polar- 
ization, the frequency of measurement was increased to a range 
of 5 - 40 MHz. Unfortunately, the attempt was only partly 
successful. As can be seen from the frequency dependence 
curves of Figure 1-23, the dielectric constant is markedly less 
dependent on frequency, and the values of permittivity are muc/i 
more reasonable than for the frequency range of 20 - 2000 kHz, 
as shown in Figure 1-19. The addition of salt, however, still 
has cin effect on the results, albeit a small effect. This 
gives one an opportunity to assess a prediction made in Chapter 
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4, where it was claimed that a one decade decrease in R 2 
should result in the C vs. frequency curve being shifted 
one-and-a-half decades to the right. R 2 » the "sample re- 
sistance," is inversely proportional to 0, the sample con- 
ductivity, while C is proportional to the dielectric constant. 

Hence, an increase of one decade in a should result in 
a one-and-one-half decade shift in Table 1-7 shows the ■ 

correspondence between this prediction and the data collected. 
The agreement is reasonable. A number of other observations 
could be made in this section concerning the permittivity, 
but it seems best to defer those remarks until the next sec- 
tion where electrode polarization does not confuse the issue, 

, Measurements in the 250 - 450 MHz Range 

The data discussed in this section were collected ac- 
cording to the procedure outlined in Chapter 5 and are tab- 
ulated in the third section of Appendix I. 

In this frequency range electrode polarization has a neg- 
ligible effect. This can be seen from the fact that the addi- 
tion of salt does not affect the permittivity and that the di- 
electric constant is not strongly dependent on frequency. 

Considering the frequency dependence in more detail, one 
notes that the changes in log e with respect to frequency, 

A Log e/Af, tend to be higher for soils higher in clay content 
but, for any given soil, A Log e/Af is independent of moisture 
content. Likewise, the change in conductivity increases for 
increasing clay content, but Aa/Af is approximately independent 
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Table 1-7. Comparison of Predicted with Actual Frequency 
Shift of C vs. Frequency Curve Due to Change 
in Conductivity. 


Predicted Frequency Actual Frequency 
Samples Compared Shift in Decades Shift'' in Decades 


9.9% 

M.C., 07o Salt 



10.0% 

M.C, 0.47oSalt 

0.606 

0.66 

12.4% 

M.C, 07o Salt 



l2.37o 

M.C., 0.4% Salt 

0.925 

1 .00 

l5.87o 

M.C, 07o Salt 



l5.97o 

M.C, 0.4 7o Salt 

0.72 

0.87 

I9.07o 

M.C, 0.27o Salt 



l8.97o 

M.C, 0.47o Salt 

Ml 

1.08 


Note: The above data are for Miami Silt Loam in the 5-40 

MHz range. M.C. refers to moisture content. 
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of moisture content. The conductivity, however, is increasing 
faster with respect to frequency in this rcinge than in the 
5-40 Mhz range. 

Now let us turn our attention tcrward the behavior of the 
dielectric constant as a function of moisture content. The 
most striking feature of Figures 29 to 34 is the linearity of 
the graphs over most of their range. This is strange since a 
simple mixing rule, for example, (fraction water) 

+ (fraction soil) , would lead one to expect the curves 

to be straight lines on a linear scale, not a logarithmic one. 
Also, note from Figure 1-34, that the curves for different 
soils have different intercepts (projected) on the abscissa 
axis. The significance of this observation can be better ap- 
preciated by looking at Table 1-8 where the permittivity of 
each soil at its wilting point is tabulated. It is possible 
to conclude from this table and Figure 1-34, that the inter- 
cepts are correlated with the wilting point moisture, implying 
that the permittivity is a logarithmic function of the avail- 
able moisture content. The logarithm of relative permittivity 
vs. available moisture content is plotted for several soils 
in Figure 1-36, bearing out the above conclusion. 

To explain these observations one can begin with the con- 
cept of availcdDle moisture. In Chapter 1 availcible moisture 
was defined as the amount of moisture in excess of the wilting 
point moisture. The availability of water to plants and the 
permittivity of water are both related to the same physical 
phenomenon, the adsorption of water. Water becomes "bound" 


Table 1-8. Permittivity of soils at the wilting point. 
Soil Relqtive Permittivity 


Chelsea (sand) 

4.05 

Crider (cl(^) 

5.07 

Crosby (silt loam) 

4.91 

Miami (silt loam 

4.^ 


f=450 MHz 
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This figure has been adapted froo Figure 1-34 iby adjusting 

the abscissa to show available moisture content rather than 

total moisture content. 
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to the surface of the soil particles and this binding makes 
it difficult for the roots to extract moisture. The rotation 
of water dipoles by an alternating electric field is also 
hindered by the binding of the water to the surface. Since 
the high dielectric constant of water is due to the ability 
of its dipoles to rotate in response to an applied electric 
field, the dielectric constant of "bound” water is much less 
than that of free water. As a first approximation, soil mois- 
ture can be thought of as consisting of two kinds of water, 
free (with a relative permittivity of about 80} and bound 
(with a relative permittivity of 5 to 10). The data from this 
investigation are consistent with this explanation, because 
moisture below the wilting point seems to have much less ef- 
fect on the permittivity. 

Examining the data in more detail, one can arrive at 
further conclusions. For exeunple, one may ask why the log- 
arithm of permittivity is proportional to the available mois- 
ture content instead of the permittivity itself being pro- 
portional to the moisture content. A possible answer is that 
all the available water is not equally availcible. Assuming 
a distribution of permittivities corresponding to a distribution 
of binding forces, one can obtain an approximately logarithmic 
relationship as illustrated in Figure 1-37. 

The top illustration in Figure 1-37 shows a hypothetical 
graph of the dielectric constant of water vs. distance from 
the surface of a soil particle. This can be used to obtain a 
hypothetical permittivity vs. moisture curve. Consider a 
reasonable simple mixing rule. 
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Figure 1-37. Effect of Transition Region. 



Distance from port icte surface 



Moisture content 



Notes The daehed curvea are hypothetical curves drawn under the 
aaaunption that the water present is either cos^letely bound to the 
soil particles (c ■ 5:^10) , or cosnpletely free (e > 80). The solid 
curves are iiypothetical curves drawn under the assumption that over a 
saiall range of distances fron the soil particles the water molecules 
are only partially hound and the relative permittivity of this water 
gradually increases from about 10 to 80 as the distance from the soil 
particle Increases. 
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^mixture'^^r soil) 

Because (fraction soil) = 1 - (fraction water), 

E =(e . -e .,) (fraction water) +e 

r water soil soil. 

Hence, the slope of the relative permittivity vs. moisture 
content curve is simply the difference between the permit- 
tivity of the water and that of the soil. By using this re- 
lationship, one can construct the relative permittivity (e^) 
vs. moisture content curve graphically as shown in the second 
illustration of Figure 1-37. Finally, the third graph shows 
the permittivity vs. moisture content curve using a logarith- 
mic permittivity scale. Superimposed on each graph is a curve 
(dotted line) showing the same quantities under the assumption 
that only free and bound water exists with no transition region. 
The figure shows that the addition of the transition region 
gives curves which are similar to the experimental curves. 

At this point, the reader is reminded that the models set 
forth here are merely presented as reasonable explanations to 
be used as a means of generating new predictions and, possibly, 
stimulating new research. In this spirit, let us consider the 
frequency dependence of the permittivity. The frequency de- 
pendence for ordinary free water is well Itnown. It is constant 
up to the "relaxation frequency" where it rapidly declines to 
a relatively small value. This relaxation frequency is above 
any frequencies used in the present work. One might hypothe- 
size that for bound water the relaxation frequency is simply 
shifted to a lower frequency. This is a reasonable assumption 

I 

since ice behaves this way. If this is the case, then one 
might conclude that the water in the transition region is 
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undergoing relaxation in our frequency range of nieasure- 

, I 

ment. This suggests a slight decline in the dielectric 
constant with frequency which should be more noticeable for 
small moisture contents. The decline should also be a little 
greater for clays, because in both cases a greater percentage 
of the total moisture is in the transition range. This is 
borne out by Figures 1-10 and 1-12. 

The conductivity data can also be interpreted in this 
s^une light. The conductivity can be considered as being com- 
posed of two components: one due to ionic conduction, the 

J 

other due to dielectric loss. As is well known, the theory 
of dipole relaxation shows that dielectric loss reaches a 
peak at the relaxation frequency. For ordinary free water, 
the dielectric loss begins to increase around 500 MHz. Since 
water in the transition region is nearing its relaxation fre- 
quency, it should show increasing loss at a somewhat lower 
frequency than does free water. This should be reflected in 
cin increasing Aa/Af , with increasing frequency which, in fact, 
does increase, as shown by the data in the Appendices. Of 
course, the conductivity increases for increasing moisture, 
because this favors increased ionic conduction. ! 

Our main interest in performing these measurements is to 
arrive at an assessment of the usefulnes of permittivity as an 
indicator of soil moisture. From these measurements one can 
conclude that the permittivity is a good indicator of soil 
moisture provided that some information about the soil is 
available. Even in the absence of any knowledge about the soil. 
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Figure 1-36 shows that the relative permittivity gives a use- 
ful estimate of the available moisture content. For the soil 
high in clay content (Crider) the correlation is least sat- 
isfactory. 

It is also of interest to note that the laboratory samples 
were very similar to the Miami samples taken from the field, 
differing only in conductivity. It is not surprising that the 
conductivities of the field samples were greater, since pore 
structure, which is absent in the laboratory prepared samples, 
increases ionic conduction. 

As a final point, let us consider the larger than ex- 
pected permittivities for small availedale moisture contents in 
clay. There are at least two possible explanations. The first 
possibility is that the trcmsition region is constructed in 
such a way that the permittivity at the wilting point is higher 
for clays. The second possibility is that soil particle sur- 
face admittance (discussed in Chapter 2) still has some small 
effect even at our highest frequencies. According to Schwartz* 
theory, (14), spherical particles 0.02 micrometer in diameter 
could have some effect up to 500 MHz. For needle-shaped par- 
ticles the effect could be enhanced. In either case it would 
seem that performing the measurements at even higher frequencies 
than 450 MHz would give better results. A higher frequency 
would lessen the effect of surface admittance and also would 
"relax" more of the dipoles. It seems possible, therefore, 
that one could find a range of frequencies for which a better 
correlation exists between available moisture content and per- 
mittivity. 
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CHAPTER 7 
CONCLUSIONS 

The results of a number of previous investigations on 
the relationship between the electrical permittivity and the 
moisture content of soils have been accounted for herein, by 
the theory of electrode polarization. Hence, the negative 
conclusions reached about dielectric methods, based on these 
earlier studies, can be disregarded. 

When the effects of electrode polarization can be neg- 
lected, it seems that an acceptable relationship exists be- 
tween available moisture content and electrical permittivity. 

A model has been developed which is based on concepts origin- 
ated by soil scientists. The model appears to explain the 
main features of the data and is consistent with the observa- 
tions. In addition, the model predicts an even better correla 
tion between available moisture content and permittivity at 
still higher frequencies than those used in the present work. 
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PART II 

Electrical Methods for Detennining 
Soil Permittivity Profiles 
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-CHAPTER 1 

STATEMENT OF THE PROBLEM 


The ultimate goal of this investigation is to explore 
the feasibility of electrical methods for determining the 
soil moisture profile, defined as the soil moisture content 
as a function of depth in the soil., In the previous part of 
this report, a satisfactory relationship was obtained be- 
tween the available moisture content and the electrical per- 
mittivity for several different soils. Thus the problem re- 
duced to determining the permittivity as a function of depth 
in the soil. 

One possible method for measuring the permittivity pro- 
file is similar to conventional techniques of determining the 
soil moisture profile. A hole could be dug or augered and an 
electrical probe inserted to measure the permittivity at 
various depths. This might offer certain advantages in cost 
and convenience over conventional methods, but electrical 
methods hold out the promise of still greater advantages. 
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Since electromagnetic radiation does penetrate the earth 
and also is partially reflected by changes in the electrical 
propertifefe of the soil, it seems possible that information on 
the permittivity could be recovered without actual physical 
contact with the soil below the surface. If a method could 
be found for accomplishing this it would offer significant 
advantages over presently available schemes in speed and 
convenience. 

In this part of the report, various remote sensing tech- 
niques for determining soil moisture content are explored. They 
are termed remote sensing techniques because actual contact 
with the object being measured is not required. To begin, 
some factors which influence the choice of methods are ident- 
ified. 

Factors Which Influence the Choice of Methods 

Electrical methods rest on the assumption that the elec- 
tromagnetic radiation penetrates the soil to a sufficient depth. 
In addition to this, particular techniques make assumptions 
about the electrical properties of the soil and other aspects 
of the problem. To discuss these assumptions and requirements 
more intelligently, one should consider the various factors 
which influence the choice of methods: electrical properties 

of the soil, distribution of the permittivity in the soil, 
temperature, penetration of electromagnetic radiation, reso- 
lution of measurements, and the kinds of data that can be 
collected. 
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Magnetic Permeability 

Some ores have high electrical permeabilities, but 
study by Lukshin et.al. (1) revealed that all the common 
soils tested had relative permeabilities of about one. Hence 
permeability can be assumed to be equal to its vacuum value 
everywhere, independent of all variables. 

For many years the Electrical Conductivity of the soil has 
been studied to determine its suitability as an indicator of soil 
moisture. The studies have established (see (26) , for example) 
that conductivity is not a reliable indicator because the ions 
in the soil strongly influence the conductivity. The quantity 
of ions in the soil vary from place to place and time to time. 

Frequently, the conductivity is sufficiently large that the 
conduction current isn't negligible compared to the displacement 
current. Although the conductivity of soil varies with fre- 
quency, it is approximately bounded by zero below and one 
mho/meter above. Actually, most soils have conductivities 
of less than one- tenth mho/meter. 

n Electrical Permittivity 

Permittivity is related to the available moisture content 
of the soil. It appears to decline slightly with increasing 
frequency between 1 MHz and 1 GHz. 

In Chapter II of Part 1 studies were cited which demon- 
strated that clays have a very high dielectric constant due to 
their "surface admittance." In the frequency range where sur- 
face admittance is important, permittivity is not a relieible 
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indicator of soil moisture. Surface admittance is the dom- 
inant effect below 10 kHz. Above 2 GFiz the permittivity de- 
clines due to the inability of the water dipoles to reorient 
themselves in response to the rapidly oscillating field. Ap- 
proximate limits on the relative permittivity of soil are 
three and fifty. 

Temperature also has some effect bn permittivity. This 
is worth noting since temperature can vary with depth in- 
dependent of moisture content. Thus permittivity variations 
could be attributed erroneously to moisture content variations 

Studies by Schofield (27) revealed maximum variations of 
10®F between soil, temperatures to a depth of four feet. Other 
studies have found differences up to 18®F. Using the formula 
given in Part I, Chapter 2, one can show that this leads to 
relative permittivity variations of five percent. If this 
error is unacceptable it should be possible to correct the 
data by using information about the temperature profile. 

Distribution of Electrical Permittivity in the Soil 

In general, the soil cannot be assumed to consist of 
layers of constant permittivity. The moisture content is 
usually a continuous function of depth, although a sand lens 
or some similar structure occasionally gives rise to a rapid 
change in moisture content. 

For most agricultural soils it is reasonable to assume 
that the permittivity, profile does not change much in the 
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horizontal direction, but there are exceptions. The pro- 
perties of alluvial soils and soils with depressions frequently 
change rapidly in short distances along the surface, but such 
areas could be avoided if desired, since they are easily i- 
dentified. 

On the other hand, it is unrealistic to assume complete 
horizontal stratification. If horizontal stratification is 
implicit in a model, then the method must be checked to in- 

I 

sure that satisfactory results can be obtained by using in- 
formation from a limited area, as is done later in this 
chapter. 

Penetration of Electromagnetic Radiation 

The electrical properties of the soil govern the depth 
of penetration of electromagnetic radiation. One measure of 
penetration is skin depth, the depth at which the electric field 
is 36.8% of its value at the surface. The usual simplified 
formula for calculating skin depth cannot be used for soil 
since the displacement current density is not negligible com- 
pared to the conduction current density in this case, or vice 
versa. A correct formula for skin depth can be derived from the 
Helmholtz equation, 

2-1 - jiovi (jwe + a) E = 0. 
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For a plane wave traveling in the +z direction. 


2-2 E = Ae 


- j /- jup (jue + o) z 


The skin depth, 6, is just the reciprocal of the real part 
of the ejponent, 


2-3 6 = 


REAl [j /- jojp (juE + a) ] 


Algebra yields 


2-4 S = 


where 


/wy /-loe + / (ue) ‘ + a" 


u = angular frequency in radians per second 
u = u u 

U^= pejnneability of vacuum (4it x 10 ^ henry/meter) . 

ii^ is the relative permeability, which is dimensionless. 

" = ^r "o 

“9 

permittivity of vacuum (1/36 tt x 10 farad/meter). 
is the relative permittivity, which is dimensionless, 
cr = conductivity in mhos/meter. 

For equal to thirty-six and a equal to two-tenths of a 
mho/meter, the skin depth is 1.12 meters at 1 MHz, 0.35 meter 
at 10 MHz, 0.2 meter at 100 MHz, and 0.16 meter at 1 GHz. 


Resolution ' 

Resolution refers to the ability of a system to distinguish 
objects wiich are close together. In the context of the 
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present discussion, the objects to be resolved might be layers 
of high permittivity separated by a layer of lower permittivity. 

A measure of the resolution of a system is the least thickness 
of the intermediate layer for which the system can distinguish 
the Lwo high permittivity layers. For a thinner intermediate 
layer, the two high permittivity layers would appear to be a 
single layer. 

Even if the soil is not layered, resolution has an effect. 
Low resolution systems tend to smooth out abrupt transitions 
and blur fine detail. In some cases this is a disadvantage, 
but occasionally the lack of resolution eliminates trouble- 
some, and unimportant , effects of thin layers and inhomogeneities 

Generally, the resolution of systems is governed by the 
wavelength of the radiation. The limit of resolution is usually 
considered to be about half the wavelength. In soils, the wave- 
length is about 1/2 to 1/7 the wavelength in vacuum. At 300 
MHz these assumptions lead to a resolution of about 8 to 25 
centimeters. The resolution of each particular system is 
discussed in the section on that system since the special char- 
acteristics of each system govern its resolution. 

Information Available 

The data for the moisture measurement system must con- 
sist of measurements made of the electromagnetic fields at or 
above the surface of the earth. This constraint rules out, for 
example, methods which involve transmittance measurements since 
there is no practical way to make these in situ. 
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General Approach to the Problem 
Assumptions 

To ifacilltate analysis it is assumed that the area to be 
measured is flat and infinite in extent. A coordinate system 
is defined as shown in Figure 2-1. The permeability is as- 
sumed to be equal everywhere to its value in a vacuum. The 
permittivity may vary arbitrarily with depth, but it is as- 
sumed constant on any plane parallel to the plane z s 0. The 
same assumption is made with regard to conductivity. The dis- 
placement current density is not. assumed to be much greater 
than the conduction current density . Sufficient penetration is 
assumed and the effect of temperature is neglected. 

The Problem ^ 

Simply stated, the problem is to determine the permit- 
tivity profile, given the electromagnetic fields at or above 
the surface. One way to approach the problem is to relate the 
known surface field to the unknown permittivity profile by 
means of Maxwell^ s equations. This involves three steps. 
First, Maxwell's equations can be written for this particular 
situation and a "formal" solution can be obtained. Since the 
permittivity is unknown the "formal" solution will contain an 
unknown function. But the formal solution will contain use- 
ful information about the electromagnetic fields in the earth. 
Next, the formal solution must be related to something that 
can be measured: data taken at or above the surface. Finally, 
a relationship must be established between the formal solution 
and the permittivity variations. 
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Within the framework of this general approach there are 
a number of techniques to calculate the permittivity profile. 
The presently available methods for measuring the electrical 
permittivity of the earth as a function of depth, have been 
considered with considerable care by one of the present authors 
(F.V.S.) and the results are contained in unpublished form. 

The two methods which seem to be best adapted to the problem 
at hand are Becher's method and Slichter's method. In the 
next two chapters these methods are considered. 
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CHAPTER 2 

MARCHENKO-SHARPE-BECHER METHOD 

This method seemed to be especially well suited for 
determining the permittivity profile of the soil, since it 
was developed for geoelectric exploration. The technique 
uses reflected electromagnetic waves as data. In this chap- 
ter, the Marchenko-Sharpe-Becher method is presented and 
its suitability for determining soil moisture profiles is as- 
sessed. Actually, it turns out that the method is not ap- 
plicable to the present problem, although at first it appears 
to be, so it is only outlined here. The mathematical pre- 
sentation is sketchy; only enough detail is being included to 
make possible an understanding of the apparent inadequacy 
of the method. 


Presentation of the Method 

The foundation of this method is Marchenko's work on 
the inverse scattering problem (28), Sharpe (29) adapted 
Marchenko's results to the problem of nonuniform transmission 
line synthesis. Sharpe also showed that a particular repre- 
sentation of the data greatly simplified the process of 


109 


obtaining a solution. Finally, Becher (30) used an analogy 
between plane wave propagation theory and transmission line 
theory to ^pply Sharpe's results to the geoelectric ex- 
ploration problem. Some extensions of Sharpe's work were 
necessary to accomplish this. 

Marchenko ' s Work 

Marchenko's analysis deals with plane waves which orig- 
inate , at infinity, impinge upon a scattering potential 
V(x), and are reflected back to infinity. For the case of 
electromagnetic waves, one can easily see that Maxell’s 
equations, applied to this situation, lead to a wave equation. 
Marchenko was concerned with quantum mechanical scattering, 
so he dealt witn the Schrodinger equation, also a wave equation. 
Since Marchenko's work has general application, the equations 
presented will not be specialized to electromagnetic terms, but 
left in the form used by Marchenko, except to express the 
Phaser as rather than e 

Although Marchenko dealt with quantum scattering, his 
development follows the general approach outlined in Chapter i 
of Part II. The basic equation is 

d^Y 2 

2-5 ^ + X y = V(x)Y, 

dx^ 

where 

Y is a function of x (e.g. electrical field intensity), 

X is a separation constant, 

V(x) is an unknown potential, and 
X is a distance. ' 
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The objective of the analysis is to determine V(x) from know- 
ledge of the asymptotic properties of Y. For convenience, 
Equation 2-5 with the boundary condition, 

2-6 Y(x,X) = 1 as x -► « 

can be transformed into an integral equation, 

2-7 Y(x,X) = e^^^ + v(t) Y(t,X) dt. 

X 

This is a type of Volterra equation, and the form of a parti- 
cular solution is well known. If one postulates that 

2-8 /"t|V(t) |dt < « 

X 

it can be shown, using a theorem of Titchmarsh, that one can 

I ' 

obtain a particular solution in the form 
2-9 y (x,X) = e^^^ +/“k(x, tle^^*^dt, 

^ X 

where K(x,t) is called the kernel. 

The relationship between K(x,t) and V(t) is discussed later in 
the chapter. The first step of the general approach outlined 
in Chapter 1 is complete; a formal solution has been obtained. 
Now it is necessary to relate this formal solution to something 
which can be measured. To begin, note that a particular so- 
lution, Yp, was found above. The general solution, G(x,X), is 
a superposition of the waves propagating in opposite directions 

G(x,X) = Y (x,X)C + Y (x,-X)d. 

r Ir 


2-10 


Ill 


This can be expressed as 

2-11 G(x,X) = ^ [Yp(x,-X) - Yp(x,X> S(-X)] Yp*(-X), 

where * denotes complex conjugation, and 

S(X) = Y^*(-X) [Y^*(X)]'>. 

P P 

S(X) is called the scattering datum and has the property that 
1 - S<X) is a Fourier transform of a function, F^Ct) , 

2-12 1 - S(X) = /“F^(t)e“^^^dt. 

starting with this information, Marchenko was able to re- 
late the kernel, K(x,t), to the scattering datum, S(X), by 
means of Parseval's equation. The relationship obtained is 
termed the fundamental equation, 

2-13 F (x+y) + K{x,y) + /”K(x,t) F (t+y)dt = 0, 0 < x < y. 

X ® 

Now the formal solution must be related to the potential, 
V"(x) . This can be done by substituting the formal solution. 
Equation 2-9, into the integral equation. Equation 2-7. 

2-14 /"K(x,t)ei^‘dt = /" V(s) 

X X X 

+ /"V(s)ds /" K(s,u)du. 

X s X 

Denote the first term of the right hcurtd side by A, the second 
by B, Using trigonometric relationships, one c^m show that 

sin X (s-x) ^jXs ^ 1^2s-x ^jXt^^^ 

sin X (s-x) ^jXu _ ju+s-x 
X x+u-s 
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Substituting these relations into integrals A and B and in- 
terchanging the order of integration, one obtains 


2-16 


2-17 


/ OO 

X 

/ oo 


jXt 


i ! V(s)dsldt, 
^ x+t J 




1 

2 


X t+x-s 


K (s ,u) du 


4 


r 

x+t 

2 


t+s-x 

V(s)ds/ K(s,u)du 

* ' 
s 


dt . 


Substituting the expression into Equation 2-14 and using the 
uniqueness of the Fourier integral representation, 

1 r® 1 t+s-x 

2-18 K(x,t) = i / V(sjds + i V(s)ds / K(s,u)du 

x+t X t+x-s 

2 

1 » t+s-x 

+ i / V(s)ds / K(s,u)du, 0 < X < t. 

x+t s 


Now let t = X and Equation 2-18 reduces to the following equation, 


2-19 K{x,x) = 7 /" V(s)ds. 

X 

To summarize Marchenko's work, let us recall that he has 

provided a relationship, Equations 2-12 and 2-13, between part 
of the formal solution, K(x,t), and a known quantity, S(X), and 

a relationship. Equation 2-19, between the kernel, K(x,x), and 

the desired quantity, V(x). Note, however, that in order to 

calculate V(x), one must first solve an integral equation, 

Equation 2-13. 


Sharpe's Work 

Sharpe dealt with the problem of constucting the charac- 
teristic impedance of a transmission line in such a way as to 
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produce a specified input admittance. By normalizing certain 
variables in the transmission line equations/ he showed that 
these equations can be put in the form of a Schrodinger equation. 
Thus, Sharpe's problem is identical to Marchenko's from a math- 
ematical viewpoint. Marchenko's results, therefore, apply to 
Sharpe's problem also. Recall that Marchenko's solution re- 
quired solving an integral equation. Sharpe sought to over- 
come this difficulty by representing the specified input ad- 
mittance as a rational function of frequency. By using contour 
integration, he was able to reduce the integral equation to an 
algebraic equation. The details are outlined below. 

Sharpe began with the transmission line equations for a 
lossless line; 

2-20 II = juL{2)I(2) , 

2-21 II = jojC( 2 )V(z) , 

where 

V is the voltage across the line, 

I is the current, in the line conductors, 

L(z) is the inductance per unit length of the line, 

C(z) is the capacitance per unit length, 

(0 is the angular frequency of the voltage and current, 
z is the distance from the input of the line. 

Define a local characteristic impedence, 2^(z), and a local 
phase coefficient, 3(z), 

2-22 Z^<z) = /lTzTTcIzT, 




I 
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2-23 6(z) = /lTzTcTTT. 

Now normalize the distance, voltage, and current, respectively. 
2-24 X = 1 (t)dt 


2-25 u(z) = V/ /Z^TzT 


2-26 v(z) = I /Z“Cz) 

Substituting these normalized functions into Equations 2-20 
and 2-21 gives , ‘ 

2-27 + p{x)u(x) - jwv(x) =0, 


2-28 


dx 

dv (x) 
dx 


- p{x)v(x) - jwu(x) = 0, 


where p(x) = j ^ In ZqZ(x). 

Eliminating v(x) in Equation 2-27 and u(x) in Equation 2-28, 
one obtains the Schrodinger equations. 


2-29 

d u(x) . 


- P (x) ] 

u(x) 


dx^ 




2-30 

d^v(x) . 
5 — 

. 

- Q(X)] 

v(x) 


dx^ 




where 

P(x) = p^ 

(X) - 

dp(x} 
dx ' 

Q(x) 


Note that P(x) and Q(x) are subject to the same restrictions 
as Marchenko's potential, V(x}. See Equation 2-8. 

Sharpe thus established the connection between transmission 
line theory and Marchenko's work. To see this, note the analogy 
between Equation 2-29 or 2-30 and Marchenko’s equation (2-5). 
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Since Marchenko's results apply to this case, one can write 
the solution to Equation 2-9 by analogy to Equation 2-9, 

2-31 u(x,u) = e^“* + /"A(x,t)e^‘^^dt, 

X 

where 

A(x,t) is the kernel, analogous to K(x,t) in Equation 2-9. 

The next step is to determine A{x,t) in terms of known quanti- 
ties. Under the assumptions that he imposed, Sharpe showed 
that the problem of determining the kernel from known quanti- 
ties , was reduced to that of solving a linear system of al- 
gebraic equations. 

Marchenko's fundamental equation, 2-13/’ is the key equation 
in the process of determining the kernel. To reduce this in- 
tegral equation to an algebraic system, Sharpe first related 
the input admittance to Marchenko's datum, which is analogous 
to a reflection coefficient. The normalized admittance of the 
transmission line is 


2-32 y(x,w) 


v{x,u) 
u (x,to) * 


The normalized admittance is 


2-33 Y(u)) = y(0,ui) . 

Under the condition that the line becomes sufficiently 
uniform as x approaches infinity, outgoing wave solutions which 
are asymptotically exponential exist and Y(l) = v( 0, A)/u( 0,1) 
is a positive real function. X is the complex frequency, 

X = w - j E. 
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It is possible to prove that u(0,u) and v(0,<i>) are 

rational functions of ui with only simple poles and that, 

knowing any one of the quantities Y(X), u(A), or v(X) 

allows one to determine the other two uniquely. So, knowing 

Y(A) allows one to determine u(0,(b). Once u(0,u) is known, 

Marchenko's theory can be used to obtain F. (t) and A<x,t) 

s 

from Equations 2-11, 2-12, and 2-13. Note that u(0, -u) /u(0 ,(ij) 
is analogous to S(X) in Equation 2-12. Now perform the in- 
verse Fourier transform on both sides of Equation 2-12. Then 


2-34 


r, _ 1 r» , u(0,-w) 

^s^^^ ” ^ ^ u(0;o)) ® 


jut 


du. 


In order to obtain the kernel, A (x,y) , one can then solve an 
equation analogous to Marchenko's fundamental equation, 2-13. 


2-35 F (x+y) + A(x,y) + /*A(x,t) F (t+y)dt = 0. 

s ^ s 

Consider this calculation in more detail. Recall that 
u(0,u) is a rational function, with only simple poles, and so 


u(0,-u) 
u (b ,u) 

2-36 


can be presented as 

u(0,-u) _ r 
" u(0,u) U-K^ 


n 


+ I , 


wheie 

i 

is a complex pole whose imaginary part is greater 
than zero, 

is a complex pole whose imaginary part is less 
than zero, 

and cr^ are complex residues. 



Contour integration over the entire upper half plane gives 
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Substituting into Equation 2-35 and performing the integration, 

n . 

2-38 A(x,t) = I f (x)e3%^ ' 

v=l 

where f is found from the system 


2-39 


-P 'I 




+ f (X) + jp',e^''v* 


Thus the integral equation reduces to a linear system of al- 
gebraic equations. Once A(x,t) is determined, one can use 
Equation, 2-31 to determine u(x,ui). 

2-31 u(x,o)) = e^***^ + /*A(x,t)e^‘^^dt 

This expression ceui be substituted into Equation 2-29 
to yield P (x) . 


2-29 


d*u(x) 


+ [lu* - P(x)]u(x) = 0. 


Upon substituting Equations 2-31 and 2-38 into Equation 2-29, 
one obtains. 


2-40 r I f (x)e^^'^ ^dt + [u^ - P(x)] 

. ^ v-l 


leJ”* + /* f f (x)ei<" * “>‘dt)=0. 

X V=1 

simplifying, using standard formulae for differentiation emd 
integration, one obtains, 


2-41 P(x) “ “2 I ^fy(x)e^‘'v^ . 

P(x) is related to Z.(x), the desired quantity. To see this. 


recall that 



118 


2-42 P(x) ^ In Ze(x))* ^2 In Z.(x). 

It is possible to show that 


2-43 


Z«(x) = 


Zo(») 


det [l-R(x) 1 
det tl-R(x) ] 


wnere 


1 is the identity matrix, 

R is a matrix whose elements are R^^ 
Zq{«) = lim Zq(x), xh-». 





* %>x 

1C 

U 


t 


Becher ' s Work 


Becher applied Sharpe's results to geoelectric exploration, 
that is, the use of electromagnetic waves to determine the 
electromagnetic properties of the earth. By exploiting cm 
analogy between plane wave propagation theory and transmission 
line theory, Becher was able to put the equations for a plane 
wave in the same form as Sharpe's Equations 2-20 and 2-21. 


2-44 

2-45 


d£ 

dz 


X (z) 


dH (z) 


= j(*>lioHy(z), dEx(z)= ju)UaHy (z) , 
dz 


=-jwe (z)E^(z) , 


where 

is the electric field intensity in the x direction, 
is the magnetic field strength in the y direction, 
is the permeability of vacuum, 
e(z) is the permittivity of the earth, 

(j is the angular frequency , 

x,y, and z are defined in Figure 2-1. 
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In addition to plane waves, Becher also treated the 
case of a raan-made source (a current distribution) and de- 
rived equations for this situation as well. Becher also 
made several extensions of Sharpe's work. First, he showed 

that the lossy line (R-L line) , where cj>>(jjc, can be cast in 

i 

the form of Equations 2-20, and 2-21, which were derived for 
the lossless line (L-C line) , where Becher then pro- 

ceeded to derive a transformation between these cases. For 
the R-L line, 

(z,w) 

= -0(2)E(2,t.l) , 

where 

a(z) is the conductivity of the earth. 

If 

2-47 Mo = Ij and jw = y*, 
then 

2-48 jtupoH(z,u) =-jYL [ j YH ( z , - jy * ) 1 . 

Define 

2-49 if(z,Y) = jYH(z,-jy®), 

]^(z,y) = E (z,-jY^) 

The result is 

2_5o =.jyL ?({z,y) , 




“3^ 
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and 

2-51 =“jYC(2) 2i(z,Y)» 

where C represents a and y represents the tremsformed angular 
frequency. Equations 2-50 and 2-51 are in the form of the 
lossless transmission line equations. 

A second extension of Sharpe's work due to Becher con- 

/ 

cerns distance normalization. Sharpe dealt with lossless 
transmission lines where the speed of propagation changes only 
slightly over distances which are small compared to the wave- 
length. Hence, his normalized distance, x, is approximately 
proportional to the true distance, z. 

2-52 X = i /^S(t)dt % - /^Bdt % ^z. 

For the geoelectric problem, however, this is no longer true. 
Becher derived a differential equation relating x and z and 
proceeded to solve it. By substituting the results into 
Sharpe's equations, Becher was able to express the variation 
of permittivity as a function of true distance. 

Requirements For Using Becher *s Method 
The source of electromagnetic energy, the data collection 
process, and the medium whose permittivity variations are being 
measured, must each satisfy certain requirements. First, con- 
sider the requirements for the source. Becher considered two 
cases, a plane wave source, and a line source. In the first 
case, the plane waves were assumed to originate from natural 
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electromagnetie disturbances due to fluctuations in iono- 
spheric currents. In general, these plane waves are very 
low frequency: usually the lowest frequency is less than one 

Hz. If plane waves are not available, one must use a line 
source. In either case, the source must have a fairly broad 
bandwidth. The bandwidth depends on the expected change in 
the permittivity as a function of depth, and can be estimated 
by exploiting the fact that the relationship between the re- 
flection coefficient, T(23), and the log of the characteristic 
impedance is in the form of a Fourier transform. 

2-53 r(2e) = [lnZ(z)]dz. 

If Z(z) is considered to be the impulse response and r(2B) is 
considered to be the bandwidth, then the product of the band- 
width, B, and the impulse response, R, of is 

2-54 BR>2. 

B and R are both measured in the radius of gyration sense. For 
an example worked by Becher, the required bandwidth is three 
decades . 

The data required by the method are the surface electric 
and magnetic fields. For the case of a plane wave source, it 
is sufficient to measure the electric field intensity vector 
and magnetic field strength vector perpendicular to- it at a 
single point. Note that both vectors are in the pleine of the 
surface. When a line source is used, however, it becomes neces 
sary to measure the fields over the entire surface. In general 
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it is possible to choose a current distribution which gives 
azimuthal symmetry. In this case# one needs to measure only 
the maugnetic and electric fields along a line passing through 
the point where the permittivity profile is desired. For a 
line source, the data are the input admittemces 

2-55 Y(0,u) = /“h (y,0,o))dy//“E (y,0,u)dy 

0 y 0 * 

for all frequencies, (The coordinates are defined in Figure 
2-1.) As a practical matter, one can truncate the integration 
at some finite value of y. This point will be explored in 
more detail in Chapter 3 of Part II. 

The Marchenko-Sharpe-Becher method is not entirely general, 
because it does not apply to an arbitrary permittivity var- 
iation. Certain restrictions are placed on the medium. To be- 
gin, it is assximed that o and £ are constant on any plane 
parallel to the surface and that a and e for large z, (dis- 
placement perpendicular to the surface), approach constants 
2 

faster than 1/z . In addition, this technique requires that 
either a>>oje or (uE>>cr, and that e or o (dependent upon which is 
the variable of interest) be constant with respect to the 
frequency, that is, tnat the medium be nondispersive. 

Conclusions 

As noted in Chapter i of Part II, the data collected in ) 
Part I show that soil cannot be regarded as nondispersive ex- 
cept possibly over a small range between a few megahertz and a 
few gigahertz. In this range, however, cje and o are of com- 
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parable magnitude so that neither iije>>a nor o>>we holds. 
Furthermore, over this range the conductivity is changing 
with respect to frequency. 

These facts imply that the Marchenko-Sharpe-Becher 
method cannot be applied to the moisture profile problem. Two 

I 

basic improvements are required: the general lossy case should 

be included in the scope of the method and knowledge of the 
surface fields at only a single frequency should be necessary. 
Both of these extensions appear to involve considereible changes 
in the development. A method which incorporates these im- 
provements was developed thirty-five years before the Marchenko- 
Sharpe-Becher method, and this technique, due to Slichter, 
will be presented in the next chapter. 
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CHAPTER 3 
SLIGHTER 'S METHOD 

Slichter's method (31) uses an antenna on the surface 
of the earth to radiate waves into the ground. Information 
collected at the surface is then used to characterize the 
electrical permittivity and conductivity profile of the earth 
vertically downward. It is believed that this method can be 
applied to the problem of soil moisture measurement. Con- 
sequently, it is discussed in considerably more detail than 
was the Marchenko-Sharpe-Becher method. 

Presentation of the Method 

A circular current sheet is located on the planar inter 
face between two semi-infinite half-spaces, in which the elec 
trical permittivity and conductivity vary only with depth, as 
shown in Figure 2-2. For reasons of convenience in analysis, 
this <()-directed surface current density, I(p), is chosen such 
that 

I (p) = p/a, 0 < p < a, 

I(p) = 1/2, “ P = a, 

I(p) = 0, 


p > a. 
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where a is the radius of the disc of current, as shown in 
Figure 2-2. This current density varies sinusoidally in 
time with cuigular frequency u. The problem is to determine 
the permittivity, e (z) , eind the conductivity, a(z), from 
measurements made at the surface, z - 0. 

The first step is to formulate the problem in terms of 
Maxwell's equations and then obtain a formal solution. 
Maxwell's equations deal with the electromagnetic field: E, 

the electric field strength, and B, the magnetic induction, 
wliich are due to the current distribution on the pleme, z = 0. 
Subscripts on these variables refer to the direction of these 
vector quantities. Symmetry implies that is everywhere 
zero. Likewise, and are everywhere zero, as can be 
seen from the following Maxwell equation, in cylindrical 
coordinates (p,<t>,z}. 



where H is the magnetic field strength. Recall that H = B/y, 
where y is the permeability, and that symmetry implies that 
operating on euiy vector is zero. Since is zero, so also 
is H,. Hence 

4 > 

f 3H 3H 

2-57 (o + 3«e)E = % 

Maxwell's equations then can be written, using the customary 
phasor notation. 
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2-58 

.A 

2-59 
2-60 

2-61 7 • S » 0 

2-62 7 • eE = p_, 

c 

where is the volume cheirge density. The boundary conditions 
at z = 0 are 

2-64 

2-65 (Bp)^2 = <®p>-Z + 

where l(p) is the - directed surface current density flowing 
in the surface z =: .0. 

From this point on, the electric field strength in the (p ■ 
direction will be denoted by E. This should not cause any 
confusion since the 4> component is the only nonzero component 
of the electric field strength. 

Substituting Equations 2-59 and 2-60 into 2-58, yields 

2-66 ^ I? " — “ 

32* 3p* P ^ p* 

By separating variables, (E = R(p) Z(z)) one obtains the Bessel 
equation 


k (^)p - Ip (!) p = 

Tz “ 3wBp 

7p 
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2-67 + i g+(A2 - 1 /pMR = 0, 

dp* . P ‘ip 

where X is the separation constcmt, and the equation 

2-68 ^ - (X* + jojpa - io*pe)Z = 0. 

dz* 

Above the surface, z < 0, where e « Sg and o = 0, this last 
equation simplifies to 

2-69 ^ - (X* - — )Z = 0, 

dz* c* 

2-70 Z = D(X) exp (±/x* - k| z) , 

where kg = u/c, and c is the speed of light. 

Returning to Equation 2-67, it is clear that R = (Xp) , 
since E and B remain finite at p = 0. Let the solution to 
Equation 2-68, which vanishes at z -*■ + '», be denoted by Z^(z,X). 
How recall that the electric field, which is the solution 
to Equation 2-66, is a product solution. In order to satisfy 
the boundary conditions, this product must be integrated with 
respect to the separation constcuit, X. 

2-71 E^. = E(z) = /"j (Xp) F (X) Z (z,X)dX, z > 0 

O 1 1 1 

2-72 E = E(z) = /“j (Xp) F (X) exp(/x* - k|z)dX, z < 0, 

“ o 1 2 

where F^ and F^ are functions chosen such that the integrals 
converge and satisfy the boundary conditions. The first of 
these conditions requires that the plus sign be used in 
exp(l/x* - k|z) for z < 0, since X* > k| for X-*«. 

= E_ at z = 0 for all X, it is clear that 


Since E 
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2-73 F (X) = F (X) 2 (0,X). 

12 11 

Now the boundary condition on the magnetic induction can 
be used to determine by using Equations 2-65 and 2-59: 


<Bp)« - <V-Z ' 


2-75 


\Tz~ J +0 \dz / -0 


ju)M,I(p} . 


Substituting 2-71, 2-72, cmd 2-73 into 2-75 yields 

„ (0,X) ^ 

2-76 { ^ 5z - 2^(0, X) /X" - ky dX = juWoKp) 

Our particular choice of l(p) makes it especially easy to 
determine F^ since a well-known integral is 


Ap/a, 0<p/a<l 

2-77 a/"j (Xp) J (Xa)dX = J 1/2, p/a = 1 

o 1 2 I 

lo, p/a>l. 


By using Equations 2-55, 2-76, and 2-77, one can conclude that 


in our case 


2-78 jauyeJ 


(Xa) = F, (X) { ■ 
2 I k 


32 (0,X) \ 

— Z (0,X) /X* - k|/ , 


for all X. 

Now the solution for E^ can be written 


2-79 E_^ = ja tavo ^ 


(Xa)j^ (Xp)Z^ (z,X)dX 


3Z^(0,X) 


- /X^ - k* Z^(o^X) 


Provided that the integrals converge uniformly, B and 

' 2 

B for z>0 are given by 
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2-80 Bp (2, p) = ^ ^ 


= ay 


J (aX}J (P^)t^Zi (z,X) 

iHiTTTi— ==r— 

‘ A* - k|z^ (0,X) 




2-81 B,<z,p) = 

z 




^ + -1 
3p p 


^ XJ (aX)Jo (pX)Z (z,X) 

”5fc'«T';x') ==== 

— - k|z (0,X) 


(3X 


Sllchter (31) shows that indeed these integrals in Equations 

2-79, 2-80, and 2-81 do converge uniformly, if 

lim e(z) < E (a constant) 

2 

Z-^-oo 

and 

lim g(z) = a (a constant). 

2 

Z-»-oo 

From the physical considerations of the present problem, these 

r 

are acceptable restraints, since the electromagnetic field will 
penetrate, effectively, no more than a few meters into the earth, 
and no exceptionally large values of e or of a normally occur 
at such depths. 

A formal solution for the electromagnetic fields has been 
found, but it involves ein unknown function, Zi(z,X). This func- 


tion must be related to some quantity which can be measured. 
Fortunately, one can relate to the magnetic fields at the 

surface. To begin, write and for z = 0 in the forms 
^ p z 


2-82 Bp(0,p) = ^“j^(Xp) kp(A)dX, 

2-83 B^(0,p) = ^“jo(Xp) k^U)dX, 


where 
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3Z (0,X) /8Z {0,X) — 

2-84 kp(X) = ayoJ^(Aa) - /[ - k| Z_ 

/az (0,A) 

2-85 k^U) = -au/J^(Aa) Z^(0,X)/f — - kl Z^i 


The objective is to invert these equations for the mag- 
netic fields in order to obtain expressions for and in 
terms of the surface magnetic fields. Now the Pourier-Bessel 
theorem states that 

o» 00 

2-86 F(r) = / Jy(pr)ydy^ F (R) (yR) RdR, 

.00 1/2 

provided J -F(R)R dR is absolutely convergent. 

B is everywhere finite and vanishes at infinity at least 
, .0 1/2 

as rapidly as 1/p . Therefore, B^p dp is absolutely con- 
vergent. is everywhere finite, but vanishes at infinity 

as 1/p in the lossless case. However, loss does in fact oc- 
cur, hence an attenuation factor proportional to e must be 
included in the expression for B . So, asymptotically, B can 
be expressed as Be~^^/p, where B is a constant. Therefore, 

,00 

/ B P dp is absolutely convergent. 

o 2» 

Expressing B. and B in terms of their real and imaginary 
P 2 

parts and, substituting into the Fourier-Bessel theoreitv one 
obtains 


2-87 b’<0,p) = /“j Up)AdX/"B! (0,C)J 

M o 1 o P I 

2-88 B".(0,p) = /"j (Xp)XdX/“B" (0,C)J. (X?) 

P O I • P I 

2-89 B'(0,p) = rJo(Xp)XdX/”B' (0,OJo(^OSd5, 

z i ' 6 z 
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2-90 B"(0,p) = /“Jo(\P)^^A/“b"( 0,^)J„ (X5)CflC, 

z o or 

where B ^ B' +jB" and B = B' + jB". 

p p^’p z z-'z 

By comparison of Equations 2-87 through 2-90 with Equations 
2-82 and 2-83, one can see that, since Equations 2-87 through 
2-90 must be valid for all values of p, 


2-91 k'(X) = A/"b'(0,p)J (Ap)pdp, 

P o P I 

2-92 k"(A) = A/"b;!{0,p)J (Ap)pdp, 

p «S P 1 

2-93 k' (A) = A/“b' (0,p) Jo (Ap) pdp, 

Z o 2 

2-94 k^(A) = aJ“b^ (0,p) Jo (Ap) pdp. 


where 


'‘p = ‘‘i 




Now define the quantity 

k^(A) 

2-95 K(X,0) = ^ 

z 


and note from Equations 2-84 , 2-85, and 2-95 that 


2-96 


K(A 


3Zi (z,Ar‘ 
“ AzV(^xT 


z = 0 

Furthermore, it now is shown that K(a, 0) can be expressed 
in terms of measurable quantities, the surface fields. By 
using Equations 2-91 through 2-94 and Equation 2-9 5, one can 


write 
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/"b‘ ( 0,p)J (Ap)pdp + j/~B"(0,p)J (Ap)pdp 

0 P 1 9 . P 1 

/"BMO,p)Jo(Ap)pdp + j/“B"(0,p)Jo(Ap)pdp 

0 * o * 


Now a relationship must be fovmd between Z ^ and the 
desired quantities, e(z) and o(z>. This relationship will 
tcike the form of a MacLaurin series expansion. As was in- 
dicated in the General Approach given in Chapter 1 of Part II, 
the starting point for finding the relationship is the sub- 
stitution of Z^ into Maxwell's equations. Recall Equation 
2-68, whose solution is Z^,^for z > 0. 

2-98 Z - (A^ + juyo - w*ye)Z = 0. 

az* » 1 

For large A , 

2-99 Z^ (z,A) = c(A)e"^^. 


Therefore, we may write 


2-100 lim K{A) 
A-«» 


azi (0,A) 



Hence, for large values of A, one can write 

CO ' ’ . 

2-101 K(A) = y a„A"", (z = 0) , 
n=0 " 

where a^ are constants, emd ae s i. 
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2-103 TT ” + (X* + jojijo - = 0. 

d z 

For large X, the asymptotic expression for v is 


2-104 v(z,X) - \l a. (z)X~”, 

n=0 ^ 

where is a function of z. It is clear from Equations 

2-101 and 2-102 that 

2-105 lim V (z , X) -*-Xk (X) , 
z-*-0 

for all X, so 

2-106 lim a (z)-+-a . 
z^O " 

Now it is possible to show that, for all z, 


2-107 a„ = 1. 


The argument to prove this is as follows. From Equationi 
2-104 we see that, for all z, 

2-108 lim v(z,X) = Xoto. 

X->oo 

Now, from Equation 2-6G, with £, y, and a bounded for all Z, and 
for finite o, we have 
9 ^ Z 

2-109 i- - [X^ + f(z)]Z = 0, 

3z^ > 

where now f(z) is bounded for all z. Hence 
2-110 lim Z (z) = Ae*^^, 
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where A is an arbitrary constant and the plus sign in the ex- 
ponent must be discarded because (z) must be finite for all 
positive values of z. Then, from Equations 2-102, 2-108, and 
2-110, it can be seen that 


2-111 a, = lim 

X-*-® X 


-lim 

X-Kn 


XZ, (z,X) 



for all z. 


To find expressions for the other a^, we proceed as follows. 
Substitute the series expansion for v into the Riccati equation 
which V satisfies, Equation 2-103: 


2-112 X[ I ctj-ci + X'^ - oj^UE + jojyo = 0. 

n=0 n=0 r=0 


Algebraic simplification yields 


2-113 


n»=-l 


-n 


da 


n+2 

n+1 r 
-3— + ) a_a. 
QZ -“n r : 


r=0 


h+2-r 


■ jajigr - a*yc 


Note that the right-hand side is independent of X , hence the 

left is also. Therefore, 

, n+2 • 

da 1 

2-114 n > 0, - 35 -t , 

da 

2-115 n = -1, = 0 3 2a a^ = 0, 

2-116 n = 0, 202 = jujpo - 
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Given a^(z), a coinplex function of z, o(z) and e(z) ceui be 

calculated assuming that u and y ao:e constant and known. 

Now, assume that a(z) and c(z) are sufficiently smooth so 

that all derivatives of interest exist. Granting this, Equa-i 

tion 2-116 shows that all derivatives of interest also exist 

for a (z). So a (z) can be expanded in a MacLaurin series. 

2 2 

du 2 r d^a 

2-117 a (z) = a (0) + 3 — z + 2. 

2 2 o dz^ 

The coefficients of this series can be determined from Equation 



O 


2-114. 

Since 

2-118 

(z) -► aj^r as z 0 , 

2-119 

a (0) = a =1, 


e 0 


da (0) 

2-120 

— i- — = 2a (0) = 2a 
dz 3 3 


d^ a (0) da (0) 

2-121 -2 = 2 - 7 ^ ■= + 2 a ^ 

d^2 4 2 

d^a (0) .2 

2-122 2 = (2a (0)a + a a + 

dz® dz® “ ^ ^ ^ 


a a ) 
2 1 



2-125 e(z) = — 2 — Re[±a (z) ] . 
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Now a (z) can be expressed as a iMacLaurin series and the coef- 
ficients of this series can be related to the kernel, K(X). 

It is worth noting at this point that K(X) can be deter- 
mined from knowledge of only one component of the surface mag- 
netic field instead of both, as was done in Equation 2-95. 

To see this, recall that 


2-126 K(A) 


-3Z (0,X) 

— 

XZ^ (0,X) 


Recall also, from Equation 2-85,' that 

2-127 k, = -apo J (^a) XZ (0,X)/ Z (0,X)-Z (0,X) (X® 

Z 2 1 I 1 



Solving algebraically for I-3Z ^ (0 , X) /3z] /fxz ^ (0 , X) ] in Equation 
2-127, yields 


2-128 K(X) 


3 Z (0,X) ayoJ (Xa) 1/2 

^ (X^ - k^) ' 

XZ (o,x) OxT X 

I 2 


One can summarize the steps involved in the calculation of 
e(z} as follows: ' 

(1) Create the current density given in Equation 2-55 and 
measure the magnitude and phase angle of the resulting 
z component of the magnetic field at the surface. 

(2) Using Equations 2-93 and 2-94, calculate k (X). 

z 

(3) Using Equation 2-128, calculate K(X). 

(4) Approximate K(X) by a power series. If X ' is taken 
as the variable, any of the many techniques for 
polynomial approximation may be used. The reader is 
referred to any standard text in numerical anaylsis. 
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(5) Using the coefficients, {a^} calculated in step 4, 
calculate the coefficients of the MacLaurin series, 
2-117, using Equations 2-114 and 2-115, as shown in 
this section. 

(6) Calculate a (z) using the MacLaurin series. Equation 

2 

2-117. 

(7) Calculate c(z) using Equation 2-125. 
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Requirements for Using Slighter* s Method 

The specifications for the source of electromagnetic 
waves are very explicitly stated in Equation 2-55. The cur- 
rent distribution used in Equation 2-55, however, was chosen 
for matheraatical convenience and does not seem to represent 
a fundamental limitation. Note that the source radiates at 
a single frequency. 

The requirements on the data collection are more strin- 
gent than those on the source. One must measure the magnetic 
fields on the surface of the earth. It is shown earlier in 
this chapter, however, that it is sufficient to measure only 
the z - component of the magnetic field. It is necessary, 
though, to measure both the magnitude and the phase of this 
component v/ith respect to the source. 

i 

One may ask how practical these requirements are. First, 
consider the measurement of phase. When Slichter developed 
the method in 1933, the measurement of phase at high frequencies 
was a very difficult task. Since then, techniques have been 
developed to perform this type of measurement. For this partic- 
ular case, it might be simplest to perform the measurement in 
two steps. First, the magnitude as a function of the radial 
distance could be measured and, secondly, the magnetic wave could 
be mixed with a wave of known amplitude and phase with respect 
to the source. The phase of the magnetic field could then be 
deduced by using the known information: both magnitudes, the 

magnitude of the sum,^and the phase of the reference. 
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Another requirement of Slichter's method is that the data 
on the magnetic field must be collected on the surface of the 
earth along a line extending from the point where the permit- 
tivity profile is desired, radially outward to infinity. These 
data are used in Equation 2-128 in the process of calculating 
the kernel , which require the calculation 

2-129 k = X/“’b (0,p) Jo (Xp) pdp. 

Obviously this requirement of making measurements of 

to infinity, is impossible. As was noted in Chapter 1 of Part 

II, the earth is seldom homogeneous in directions parallel to 

the surface of the earth for long distances. It is also clear, 

however, that the integral in Equation 2-129 converges rapidly, 

since both B and Jo(Xp) decrease rapidly with increasing p. 

k should be calculable to a good approximation by using data 
z 

collected within only a circle of radius of the point where 
the moisture profile is desired; 

2-130 k^ - (0,p) Jo (Xp) pdp. 

In order to obtain an indication of the required magnitude 

of an example was worked. The procedure consisted of; 

choosing a permittivity profile, calculating the exact kernel, 

using this to calculate the resulting exact surface magnetic . 

field, then attempting to recover the kernel by using Equation 

2-130 for various values of p . Figures 2-3 to 2-6 show the 

m 

results, which are discussed later. 

Consider the steps of the calculation in more detail. 
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Figure 2-3. Real part of exact magnetic field. (B ') as a function 
of distance from the center of the antenna. 
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Actual values 

Calculated values using Surface Reid 
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Figure 2-6. Imaginary part of J2(.1A) Z^CO.X) / aZifO.X) -Zi (O.A) A”^-kg 
versus the separation constant, A. This is the imaginary part of k as given 
by Equation 2-85, divided by (-apoA). 



145 


(1) A constant electrical permittivity profile was chosen 

to faciliate calculation. The actual value chosen 
for was 4.5. 

(2) was calculated by using the equation below 

B Z 

2-131 k, = ap;^J (Aa)Z (0,A)/t^(0,X) - /A*- k^Z (0,A)] 

Z 2 1 dZ ol 

In this particular case Z^ can be calculated easily 
from Equation 2-68, because it is being assumed that o = 0 
and that e(z) = 4.5, a constant. The function, k^ , can 
be calculated exactly for each value of A. 

(3) Using the calculated exact values, for k , the z - com- 

z 

ponent of the surface magnetic field was calculated ac- 
cording to Equation 2-83. 

2-83 B^(0,p) = /"je (Ap)k^ (A)dA. 

Z o z 

B , of course, could not be evaluated exactly, since 

Z 

the integral in Equation 2-83 is not tabulated. It will 
be seen, however, that for the particular k^ in this 

i 

excimple the integral converges rapidly. This enables 
one to estimate the quantity by calculating successive 
integrals for increasing values of the upper limit and 
checking the agreement. The difference between successive 
integrals is defined as D. 

(Ap)k^(A)dA - ^^"+^J,(Ap)k^(A)dA| 

|'/’^^Jo(Ap) k^(A)dA| 

O ^ 
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When this quantity became less than 0.005, the es- 
timate of D^(0,p) was accepted. This limit of 1/2% 
for D is consistent with the accuracies expected from 
this method of measuring soil moisture. The field 
components are sketched in Figures 2-3 and 2-4. 

(4) Then a quantity proportional to k was calculated just 
as it normally would be in a practical implementation 
of Slichter's method using Equation 2-130. 

2-130 == (0,p)Jo(Xp)pdp. 

The values of p^ used were 0.5 meter, 1.0 meter, 1.3 
meters, and 2,0 meters. The results are plotted in 
Figures 2-5 and 2-6. 

The reader may well be curious concerning the sig- 
nificance of these results. As was explained pre- 
viously, the calculation of k , and then K(A) is a key 
step in the calculation of the permittivity. So, if 
it were impossible to calculate the kernel, Slichter's 
method would not be practical. The primary physical 
limitation in calculating the kernel is the lack of 
homogeneity in any plane parallel to the surface. Hence 
it is essential to verify that the kernel can be cal- 
culated by using data from within a limited surface area. 
The results thus far indicate that tliis can be done. 

Another key step in the calculation depends upon 
expanding e(z) and o(z), the permittivity and conductivity^ 
in MacLaurin series. Since t (z) and a{z), are physical 
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quantities, it is clear that step and oscillatory 
discontinuities are not possible. While this sat- 
isfies the purely theoretical requirements for a 
MacLaurin expansion, practical constraints dictate 
stronger conditions. Specifically, e(z) and a(z) 
should be reasonably smooth functions of z, requiring 
few therms for an accurate expansion. 

As it turns out, the inherent limitations in 
resolution of about one-half wavelength, or about 
5 inches or 12 cm (a commonly accepted limit for 
resolution in instruments utilizing wave phenomena) 
help to achieve this result. Spatial variations in 
permittivity and conductivity which are significantly 
smaller than the wavelength have little effect on the 
the reflected wave and are indistinguishable at the 
surface from smooth transitions. This artificial 
smoothing, due to the finite resolution of> the system, 
causes no difficulty to soil scientists who seldom re- 
quire fine resolution. 

Review of Slichter’s Method 

The calculation of permittivity using Slichter's method 
can be outlined as follows: 

(1) Generate ^ electromagnetic field by using a magnetic 

dipole amtenna with the current, distribution: 

{ p/a, 0<p<a, 

1/2, psa, 

0, p>a. 
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The frequency should be of the order of 300 MHz , 
as shown in Part I. 

(2) Measure the magnitude and phase with respect to source, 
of the z - component of the magnetic field along any 
line extending radially away from the antenna. This 
must be done out to a distance of about one wavelength 
(approximately one meter) from the center of the antenna. 

(3) Calculate the kernel using the equation below. 


2-132 KU) = 


avi„J (Xa) 


X/'^mB (0,p)J„(Xp)pdp 


/X^- k| 


(4) Expand the kernel as a power series. 


2-133 K(X) - a X 
n=0 ^ 


(5) By exploiting the relations outlined in this chapter, 
calculate the coefficients of the ilacLaurin series ex- 
pansion of a,^. 


2-134 (0) = a , 

2 

1 (0) = 2a , 

2 3 

02 ( 0 ) = 2a + a^ , 

4 2 

etc. 

(6) Calculate a (z) . 

2 

a (0) 

2-135 = a|0) + 012 ( 0 ) z + ^ ■■ z^ 

( 7 ) Calculate the permittivity. 

2 


2-136 ^ - 


Re [a (z)]. 
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Conclusions 

Slichter's method appears to be feasible, since its as- 
sumptions and requirements seem to be consistent with the pro- 
perties of soil reported in Part I of this thesis. In addition, 
modern technology can satisfactorily perform the operations re- 
quired. 

More work is required to implement the method. It would 
be useful to calculate a number of examples using different 
permittivity profiles. In addition, one should explore the pos- 
sibility of relaxing the constraints on the current distribution 
(Equation 2-131) . Perhaps the first extension should be an in- 
vestigation of whether or not the method can be used if the cur- 
rent distribution consists of a circular ring of current. This 
would be one of the easiest current distributions to generate 
and one of the easiest about which to measure the resulting mag- 
netic field. The questions of resolution in depth, accuracy of 
the results, and actual limit of the depth to which the moisture 
content can be measured, should also be investigated. 
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APPENDIX I 
DATA TABLES 

Preliminary Investigation 

The tcUsles in this section give the measured capacitance 
resistance, permittivity, and conductivity for various sam- 
ples of soil. A complete discussion of how the data were 
obtained can be found in Chapter 3 of Part I. 


pUBCEDli^G 


pagb 


not filmfd 






Crider Clay 30.4% Moisture Content Crider Clay 31.7% Moisture Content 


156 


u 

u 

c 


m 






to 

VO 

m 

fsl 

o 


in 

o 

p 

VO 

VO 

VO 

VO 

m 

in 

in 

cQ c: 

• 

• 

• 

• 

• 

• 

• 

■H a; 
n 
0) 
os 

rH 

rH 


rH 

rH 

rH 

rH 


o 

s 


p 

■rH 

U h 

O 

m 

o 


m 

• 

VO 

o 


IQ a 

m 

r' 


VO 

OV 


in 

a 

IQ 

u 

rH 

rH 


OH 

iH 





>1 

0 

c 

V 

9 N 

era 

O 

O 

o 

o 

o 

o 

O 

01 X 

CN 

in 

o , 

o 

o 

o 

o 

p 

a 



f— 1 

CM 

m 

o 

i-H 

o 

r>H 


a 

a> 

•p 

c 

o 

u 

Q> 

9 

P 

n 

•H 

o 

s 

tfp 

00 

n 


P 

0) 

•o 

•H 

M 

u 


c 


■V 

P 

o 



M 

00 



n a 

fM 




o 

o 

ov 

■H X 

• 

• 

4 

t 

• 

• 

• 

n 

CH 

fM 


M 

fM 

M 



V 

u 

c 

IQ 

P 

•H 

U h O 

iQ Oj in 

a o 

lO r-l 

u 


o in ^ 00 p» m 

(M V in r* lO 

^ (N iH 


>1 

u 

c 

Q> 

9 N 

CTI O 
Q) ^ fM 
p 

a 


O O O O O 

in o o o o 

r-l <N in O 


0) 

u 

c 

IQ 



in 




m 

p 

c 

<u 

u 

9 

IQ 

in 




in 

m 


4J 

o 

00 

VO 


fM 

00 

fO 

0) 

P 

VO 


fM 

o 

VO 

ro 

ov 

u) c: 


VO 

VO 

VO 

VO 

in 

in 

p 

(0 c: 

VO 

VO 

VO 

VO 

m 

in 


•M M 

• 

• 

• 

t 

• 

• 

• 

a 

■H 

• 

• 

• 

• 

• 

• 

• 

u 

0) 

a 

. rH 

rH 

rH 

rH 


rH 

rH 

0 

u 

U) 

0) 

X 

rH 

rH 

rH 

rH 

rH 

rH 

rH 


p 

9 


0) 








P 

0> 







u 








u 

u 







c 








■p 

s 







IQ 








Q 

IQ 







P 






in 


X 

P 





in 

m 

•ft 






• 



•H 



1: 


• 

• 

u a 

o 

m 

m 

00 

o 

<N 

VO 

dP 

u a 

O 

o 

m 

m 

fM 

r* 

IQ a 

rH 

00 

r- 

VO 

o 


m 

VO 

iQ a 

o^ 

ov 

oo 

r» 

O 


a 

CN 


IM 

rH 

r-l 



• 

a 

rH 


CM 

iH 

t-H 


IQ 

rH 







fM 

iQ 

rH 






U 








n 

u 








>1 








rH 

fp 








u 








u 

u 








c 









c 








0) 








p 

0> 








9 N • 








V 

9 N 








O’ a 

o 

O 

o 

O 

O 

o 

o 

'O 

era 

o 

O 

o 

o 

O 

o 

o 

0) X 

fM 

m 

o 

O 

O 

o 

o 

•H 

Q> X 

CM 

in 

o 

o 

o 

o 

o 

p 



rH 

fM 

in 

o 

o 

P 

p 



rH 

fM 

in 

o 

o 

a 






rH 

fM 

u 

a 






rH 

CH 


2000 




Clay 41.3% Moisture Content Crider Clay 42.7% Moisture 


15B 


o 

o 



c 









in 

o\ 

VO 



00 



c 

■p 

m 

m 

in 

m 



n 

0) 

(0 a 

(N 

IN 

IN 

(N 

(M 

CN 

CM 

■4J 


. 

. 

. 

• 

• 

• 

. 

a 

n 

o 

o 

o 

o 

o 

o 

o 

0 

0) 








u 

Oi 









0 ) 

u 


fi 

■p 

•H b 

m 

o 

o 

a\ 

VO 

• 

o 

m 

• 

rH 

U 04 

VO 


o 

CM 

00 


VO 

m 

04 

m 

u 

00 

cn 

CM 

H 





>1 

u 

c 

O N 

3 a 

o 

o 

o 

o 

o 

o 

O 


(M 

m 

o 

o 

o 

o 

o 

0 

p 

b 



pH 

CM 

in 

o 

H 

o 

CM 




(Held in container longer than usual] 



Crider Clay 43.4% Moisture Content Crider Clay 45.4% Moisture Content 


159 


c 


<0 

o 

GO 

VO 

if) 

CN 

a\ 


4J 

CM 

rH 

rH 

H 

rH 

o 

O 

n c: 


rsi 

<N 


CN 

CM 

CN 

■ri J< 

• 

• 

• 

• 

• 

• 

• 

a 

0) 

DO 

O 

o 

O 

o 

O 


O 


0) 

u 

g 

■p 

•H b 

m 

o 


n 

n 

m 

• 

GO 

u a 

'a* 

rH 

GO 

CN 

GO 

VO 

if) 

10 

b 

Id 

u 


m 

rH 

iH 





& 

c 


0) N 

3 X 

O 

O 

o 

o 

o 

o 

O 

O'M 

CN 

in 

o 

o 

o 

o 

o 

0) 

M 

b 



rH 

(N 

in 

o 

rH 

o 

CN 


u 

c 

3 

o 



<N 

CO 

ry 


■P 

in 




cn 

n 

CN 

a X 

IN 

CN 

CN 

CN 

CN 

CN 

CN 

-H M 

• 

• 

• 

• 

• 

• 

• 

a 

O 

o 

O 

o 

o 

o 

O 


I 


0) 

u 

s 

*i 

■H b 

o 

CN 

o 

00 


00 

in 

« 

00 

U b 


CN 

a\ 

CN 

GO 

VO 

m 

3 

s* 

U 

r* 

cn 

rH 

rH 





>1 

u 

c 

4) N 

3 X 
D'Jrf 
01 

b 


o 

o 

o 

o 

O 

o 

o 

fN 

in 

o 

o 

o 

o 

O 



i-H 

!N 

m 

o 

o 






rH 

CN 



■p 

0) 





c 

u 





(1) 

c 





4J 

(8 





c 

•P 

ro 

O 

o^ 


0 

(0 c: 



m 

m 

o 

•H j<; 

• 

• 

• 

• 


u 

fH 


rH 

rH 

0) 

0) 





M 

C8 





3 






•P 






U) 






■H 

0) 





0 

u 





X 

c 






iQ 


m 



OP 

4J 


• 



tT 

•H h 

O 


m 

rn 

• 

U & 

Ln 

00 

lO 

m 

00 

(8 

iH 





a 





i 

i8 





<3 

5 

U 





M 

■P 






iH 

>1 





•H 

u 






c 






0) N 





•H 

3 S 

O 

o 

o 

o 

g 

KTX 

CM 

m 

o 

o 





rH 

fN 

■H 

M 





s 

h 









4J 

0) 




3 

u 




(U 

c 




■P 

3 




c 

•P 

o 

fH 


0 

m c: 

m 

CM 

rH 

u 

■H ^ 

• 

• 

• 


D) 

rH 

fH 

rH 

0) 

0) 




p 

DS 




3 





■P 





n 





•H 





0 

0) 




X 

u 

e 


in 


dP 

§ 


• 

• 

m 

4J 

rH 

iH 

'S* . 

a^ 

•H a 


m 

CM 

• 

u a 





C8 





a 




g 

18 




18 

u 




0 





I-) 





•p 





fH 

Sh 




•H 

o 





a 





0) N 

o 

o 

O 

•H 

3 a 

o 

o 

o 

1 

CTJ< 

m 

o 

o 

CO 

V 



CM 

•H 

M 




S 

h 


in 

69 

lO 


74 

44 

• 

• 

• 

• 

• 

• 

CM 

CM 

CM 

iH 

o 

o 


VO 

m 

in 


in 


• 

• 

• 


• 


00 

M* 

cn 

CN 

m 

■ M" 

pH 

CO 

M* 


cn 

CM 


O O O O 0-0 o 

(N in o o o o o 

rH <N in o o 

fH (N 


Miami Silt Loam 10.1% Moisture Content Miami Silt Loam 10.8% Moisture ConteAt 


161 


•u 

in c: 

o\ 

at 

as 

VO 

iH 

in 

75 

•H 

• 

• 

• 

• 

• 

• 

• 

n 

as 

as 

as 

as 





CO 



35 

* 

« 

• 

• 

as 


m 

n 


>d 

•M 

•rl h 

o 

• 

o 

, » 

in 

• 

o 

• 

o 

• 

m 

• 

o 

' • 

o 

• 

in 

id 

■p 

•H h 

O 

o 

<N 

in 

• 

in 

in 

• 

pH 

n 

m 

U Q. 

o 

VO 

in 

r* 



pH 

pH 

u a 

(N 

o 

CN 

r* 

in 


m 

(Q 

a 

« 

u 


in 

m 


in 


fO 

0 

« 

a 

<0 

u 

r> 

<N 

pH 






3 S 

o 

O 

o 

o 

O 

o 

o 

0*Ji< 

CN 

in 

o 

o 

o 

o 

O 

V 



rH 

«N 

in- 

o 

o 

p 






pH 

CN 


V 


O O O 
O O O 
(N m o 


in in m m 









dfi 

a 






m 


in 

in 

as 

■p 






■ 


• 

• 


b4 

O 

00 

pH 

00 

o 

o 

pH 

in 

VO 


u a 

in 

ro 






rn 

CN 

pH 

<a 

CN 

pH 


iH 


o 

O 

o 

O 

o 

o 

o 

•H 

0) N 
P X 

o 

O 

o 

O 

O 

o 

CN 

m 

o 

o 

o 

o 

o 

s 

O'M 

CN 

m 

' o 

O 

o 

o 



pH 

CN 

m 

o 

pH 

o 

CN 

i3 

rH 

s: 

Q) 

u 

bt 



pH 

CN 

m 

o 

pH 



2000 27.5 5.3 2000 


silt Loam 17.0% Moisture Content Miami Silt Loam 17.1% Moisture Content 


162 


0) 

u 

% 



iH 


00 

n 


c 

u 

§ 

*i 

m 


m 

(N 

00 

IN 

rH 

0) 

■p 

(0 cz 

cn 


m 

m 

IN 

CN 

rH 

■u 

n c: 

•H Jit 

• 

• 

• 

• 

• 

• 

• 

c 

•H ,K 

m 

<i> 

P6 


rH 

pH 

rH 

rH 

rH 

rH 

0 

u 

d) 

u 

to 

0) 

a 


b 

0. 


V 

O 

c 

la 

4J 

•H 

U 

(0 

Cu 

m 

u 


>1 

u 

c 

0> N 

3 K 

CP^ 

O 

U 

h 


o 

CM 


O 


(N 







t-q 







•P 







iH 







•H 

o 






CO 

c 

(1) N 

O 

o 

O 

o 

O 

•H 

3 a 

o 

o 

o 

o 

o 

s 


rH 

IN 

m 

o 

o 


o> 




rH 

CN 

•H 

p 







a 


00 


CO 


in 

VO 


a 

VO 


m 


0) 

u 

c 







in 

in 


Id 

m 


n 

m 

in 

in 

o 






• 

. 

dP 

jj 

. 

• 

• 

• 

. 

•. 

• 

m 

iH 

m 

r- 

o 

m 

ov 

VO 

•H a 

00 


o 

00 

VO 



(N 


m 




CN 

. 

u a 

m 


fO 

CN 

CN 

CN 

CN 


Id 

a 

Id 

u 


o 

fM 


O 

irt 


o 

o 


o o o 
o o o 
(N m o 


§ 






m 


•p 

c 

3 

Id 


in 

o 

00 

CN 

rH 

4J 

GO 

in 

m 

00 

VO 



0) 

■p 

00 


C^ 


r' 

VD 

u a 

ro 

m 


CN 

o 


CN 

■p 

m a 







■H Jt 

• 

t 

• 

• 

• 

• 

• 

3 

•H 

• 

• 

• 

• 

• 

• 

U) 

a 

CN 

CN 

CN 

CN 

CN 

rH 

rH 

0 

u 

<u 

u 

0) 

a 

o 

o 

o 

o 

o 

O 


<u 

u 

c 

Id 

+j 

T» CD 

o a 

Id 

a 

Id 

a 


>1 

u 

c 

0) N 

3 

(U 

h 


tM 

in 


o 


IN 

OV 


o 

CN 


o 

in 


o 

o 


o 

o 

(N 


o 

o 

in 


o 

o 

o 


o 

o 

o 

(N 


Q> 

O 


0) 

u 


m 

n 


in 


i9 




in 

m 

in 

• 

. 


. 

dP 

■p 




• 

• 

• 

in 

00 

00 

ov 

o> 

•p a 

O 

o 

OV 

ro 

00 

rH 

CO 


VO 


• 

u a 

rH 


r* 

in 

cn 

n 





CO 

Id 


rH 






a 

Id 

o 


>1 

0 

c 

0) N 

3 a 

a> 

a 


o 

fN 


O 

in 


o 

o 


o o o 
o o o 
(N m o 


2000 28.0 0.450 2000 



80% Kaolinite/20% Chelsea Sand-33. 6% M.C> 80% Kaolinite/20% Chelsea Sand-35. 9% M.C 


163 


■u 







VO 

m 

■p 

u c: 







• 

1 

to c; 

•H 

m 

m 



rH 


fH 

-d 

•H 

to 

00 

00 



in 

ro 

f>4 


to 





r* 


in 


\o 

X n) 

p 



• 

• 

• 

• 

• 

y p 

-H h 

rn 

o 

00 



eo 

CN 

•H 

u a 

CN 

CN 

rH 



rH 

i-H 

dP 0 


VO VO in m 

fM CN "CN tN (N (N 


O 

O 

o 

o 

o 

o 

O 

rtJ 

(U N 

3 X 

CN 

m 

o 

o 

o 

o 

o 





rH 

<N 

in 

o 

rH 

o 

CN 

80% 

(U 

p 

fa 


o o 
(N in 


o o o 
o o o 
fs in o 


m c: 

CO 


m 

n 

c\ 

<N 


•H M 

• 

• 

• 

• 

• 

• 

« 

to 





VO 

VD 

in 




o 

o 

o 

O 

o 

32.9 

8 

8 

o 

rH 

00 

rH 



® H r» 

ro <s fM 


m 

in 

in 

in 

00 

• 

• 

• 

• 

• 


VO 

VO 

VO 

m 

CN 

CN 

CN 

CN 

fM 


0> N 
3 X 

o 

o 

o 

o 

o 

o 

o 

(0 


fM 

m 

o 

o 

o 

o 

o 


(U 

U 

fa 



rH 

CN 

m 

o 

rH 

o 

CN 

VO 

\ 

CO 


' C 

0) N 

3 X O 

tM 

0) 


o o o 

o o o 

fM in o 



2000 21.0 1.32 2000 


164 


o 

s 


dp 

VO 



OP 

o 

00 


0 ) 

U 

S 

+J 

n 

•H 

(0 


c: 


V 

u 

i 

•p 

■H 

o 

ID 

CX 

10 

u 


h 

a 


>1 

o 

c 

V N 

a s 

CTJC 

<i> 

p 

h 


8 8 8 8 0 0 0 
o o ^ 


o 

m 

in 

m 


o 


• 

• 

• 

• 

• 

* 

• 


VO 

VO 

VO 

VO 

VO 

in 

CN 

CM 

CM 

CM 


fN 

CN 


o 

O 

o 

O 

o 

o 

o 

(N 

in 

o 

o 

o 

o 

o 



rH 

CN 

U1 

o 

o 






r-i 

(N 



Muck Soil 9% Moisture Content Muck Soil 11.1% Moisture Content 



2000 2000 4.45 1,100 



Muck Soil 11» 7% Moisture Content Muck Soil 13.6% Moisture Content 


166 


0) 


u 

ti 

lO 

■p 

cB c: 

o 

o 

o 

o 

•H 

o 

o 

o 

o 

u 

o 

o 

00 


V 

OS 

r-i 





o 


o 

rH 


VD 


0) 

u 

s 

•H fa 
U 0. 
<0 

(0 

u 


>1 

o 

c 

V N 

9 S 

vx 

Q) 

u 

fa 



CM 


lA 


00 

ir 

•p 



m 


in 



m 

» 

• 

• 

• 

• 

• 

• 

01 



• 


• 



• 

m 

ro 

iH 

O 

o\ 

00 

00 

•H 

-H fa 




fH 

iH 

in 

o 

iH 

•H 

iH 

iH 




0 

u a 

00 

VO 

in 


m 

CM 

CM 


o 

CM 


O 

in 


o 

o 


dP 

00 


in 


o 

s 

■p 

VO 

25 

o\ 

VO 


o 

o 

0) c: 

• 

« 

• 

• 

• 

• 

a 

•H ^ 



m 

n 

O 

00 

VO 

01 

1— i 

rH 

rH 

iH 

iH 




0 

u 


10 

a 

<0 

u 


& 

a 






CO 

0) N 








o 

o 

o 

o 


9 S 

o 

O 

o 

O 

O 

o 

o 

o 

o 

o 

o 

AS 


fN 

in 

o 

o 

o 

o 

o 

N 

m 

o 

o 

O 

0) 



rH 

(N 

m 

o 

o 



crJ 

CM 

9 

P 






iH 

m 





S 

fa 









0) 

u 

c 

<0 

+» 

n 

•H 

n 

0) 

os 


c: o 

X o 

o 


o 

o 

o 


0) 

u 

% 

4J 

•H fa 

u a 

lO 

u 


>1 

u 

a 

0) N 
9 M 
VX 

a> 

u 

fa 




<Tl 

m 


o 

M 


0) 

o 








13 





m 

00 







-p 





a 

• 

o 

o 


o 



0 ) c: 

o 

o 

o 

O 


CN 

o 

o 


o 


■p 

•H X 

o 

00 

in 

o 

in 

fO 

o 

o 


in 


a 

(0 

in 

cs 

rH 

iH 



iH 

iH 

















■p 

Oi 












c 













0 













u 













0) 

V 












p 

o 












9 

q 







o 

VO 

m 

m 

iH 

•P 

(3 

o 

o 

in 

VO 



• 

• ' 

• 

t 

• 

a 

■p 

a 

a 

a 

a 

a 

• 

m 

CO 

<n 

m 

CO 

•H 

-H fa 

fS 


ro 

o 

VO 



iH 

iH 

iH 

rH 

0 

u a 

n 

nj 

CM 

OJ 

iH 

rH 






S 

Id 








VO 


o 

o 

O 

O 

o 

o 


01 N 
9 X 

m 

o 

o 

o 

o 

o 

X 



H 

CM 

m 

o 

o 

u 

(1) 





H 

CM 

9 

2 

P 

fa 


u 


& 

c 


o 

CN 


O 

m 


o 

o 




O 

o 

in 


o 

o 

o 


2000 



Muck Soil 16% Moisture Content Muck Soil 25.2% Moisture Content 


167 



2000 20.3 2.69 2000 70 



168 


I 



c 









<0 









■p 


o 

o 

o 

O 

o 



n c: 


o 

o 

o 

O 

o 

o 

■p 

•H 

8 

o 

o 

o 

o 

N" 

VO 

c 

n 




«» 



p* 

Q> 

0) 


o 

o 


N* 

iH 


■P 

Oi 


iH 

iH 





C 









0 









u 









(1) 









p 









3 

0) 









o 








m 

3 








•H 

(0 








0 

■p 

m 

o 

IN 

IN 

in 

O ' 

in 

s 

■H b 


VO 

VO 

in 

N- 


n 


u a 

• 

• 

• 

• 

• 

• 

• 

dP 

Id 



N* 

N> 




CO 

a 








• 

Id 








in 

o 








iH 








- 

i 

0 

>1 









u 









c 








■p 

01 N 








fH 

3 a 

o 

o 

O 

O 

O 

o 

O 

•H 

crx 

CM 

m 

O 

O 

O 

o 

O 

00 

(U 



r-l 

(N 

m 

o 

O 


p 






fH 

fS 







170 . 


V 



u 



in 

ON 

in 




■p 



• 

• 

• 




n c: 

CNI 

OV 

'S' 

m 





•H 

m 

'S' 


m 

CM 

(N| 



CO 









0) 








■p 

« 








c 









(U 


















c 

0) 








0 

o 








u 

c 









<0 








0) 

■p 

m 

o 

VO 

o 

iH 

n 

r> 

u 

•H h 

• 

• 

• 

• 

- • 

• 

• 

3 

u a 

m 

OV 

VO 

n 


o^ 

CO 

■P 

<a 

tN 

iH 

1 — 1 

rH 

iH 



U] 

Cu 








■H 

tt 








0 

u 








s 









dP 









VO 

>1 









o 








• 

c 








VO 

0) N 
3 K 

o 

O 

o 

o 

O 

o 

o 

'0 

CTJ^ 

(N 

in 

o 

o 

o 

o 

o 

c 

0) 



1— 1 

ns 

to 

o 

o 

nj 

P 






H 

(N 

ui 

b 










171 


Properties of Miami Silt Loam 
In the 5-40 MHz Range 

These data were collected using the General Radio 1606 
R-F Bridge. The procedures used to prepare the samples and 
collect the data are presented in Chapter 5 of Part I. A 
description of Miami Silt Loam is given in Appendix II. 

The symbols used in the tables are as follows: 
f s frequency in MHz 

x*f = the product of reactance in ohms and frequency in 
MHz 

R = resistance in ohms 
e = relative permittivity 
Log e = logarithm to base 10 of e ' 
o = conductivity in mhos/meter 
M.C. = moisture content 




Miami Silt 

Loam, 9.1% 

M.C. 


f 

x*f 

R 

e 

a 

Log e 

5 

391 

148 

9.91 

0.0055 

0.996 

10 

920 

90.0 

9.84 

0.0057 

0.993 

20 

1685 

44.7 

8.07 

0.0051 

0.907 

40 

1965 

15.7 

7.86 

0.0058 

0.895 
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Miami Silt Loam» 9.9% M.C. 


f 

x*f 

R 

e 

0 

Log e 

5 

287 

103 

14.9 

0.0077 

1.17 

10 

571 

74.0 

11.6 

0.0088 

1.06 

20 

1020 

45.0 

9.69 

0.0100 

0.98 

40 

1620 

20. 0 

8.47 

0.0096 

0.92 




Miami Silt Loeun, 

10.3% M.C. 



f 

x*f 

R 

e 

0 

Log e 

5 

556 

140 

12.4 

0.0046 

1.09 

10 

1000 

88. 0 

9.99 

0.0052 

0.99 

20 

1560 

42.0 

8.69 

0.0055 

0.93 

40 

2030 

15.3 

7.68 

0.0054 

0.88 


f 

x*f 

Micuni Silt 
R 

Loam, 10.5% M.C. 

, t 

e a 

Log c 

5 

360 

121 

13.0 0.0064 

1.11 

10 

703 

80.0 

11.0 0.0073 

1.04 

20' 

1185 

46.5 

9.16 0.0084 

0.96 

40 

1782 

18.0 

8.23 0.0077 

0.91 



173 




Miami 

Silt Loam, 

12.4% M.C. 


f 

x^f 

R 

£ 

0 

Log 

5 

404 

127 

12.8 

0.0059 

1.10 

10 

826 

82.5 

10.7 

0.0063 

1.03 

20 

1360 

42.0 

9.3 

0.0068 

0.970 

40 

1836 

13.9 

8.5 

0.0059 

0.930 


f 

x*f 

Miami Silt Loam, 
R e 

15.8% M.C. 

0 

Log £ 

5 

176 

78.0 

17.4 

0.0111 

1.24 

10 

420 

55.0 

15.8 

0.0119 

1.20 

20 

797 

29.0 

14.6 

0.0122 

1.16 

40 

1130 

10.0 

13.4 

0.0106 

1.12 




Miami Silt 

Loam, 

18% M.C. 


f 

x*f 

R 

e 

a 

Log E 

5 

442 

120 

14.3 

0.0056 

1.15 

10 

812 

69.0 

12.8 

0.0063 

1,10 

20 

1220 

31.8 

11.4 

0.0069 

1.05 

40 

1578 

10.0 

10.2 

0.0059 

1.00 


174 


Miami Silt Loam^ 21.2% M.C. 


f 

x*f 


R 

£ 

o 

Log e 

5 

143 


. 65.0 

20.7 

0.0134 

1.31 

10 

302 


48.0 

17.0 

0.0155 

1.23 

20 

600 


30.0 

14.8 

0.0170 

1.17 

40 

995 

j 

12.3 

13.8 

0.0151 

1.14, 


f 

x*f 

Miami 

R 

Silt Loam, 

£ 

24.1% M.C. 

a 

Log e 

5 

334 

95.0 

IB.O 

0.0074 

1.25 

10 

491 

56.0 

16.0 

0.0105 

1.20 

20 

475 

28.5 

13.6 

0.0091 

1.13 

40 

1460 

8.20 

11.1 

0.0057 

1.04 




Miami 

Silt Loam, 

9.5% M.C. , 0.2% 

Salt 

f 

x*f 

R 

E 

a 

Log e 

5 

130 

73.0 

15.6 

0.0127 

1.19 

10 

311 

57.6 

13.0 

0.0140 

1.11 

20 

665 

39.5 

11.0 

0.0152 

1.04 

40 

1200 

19.6 

10.0 

0.0148 

1.00 
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Miami Silt 

Loam, 10. 3% 

M.C., 0.2% Salt 


f 

x*£ 

R 

e 

CT 

Log e 

5 

261 

108 

13.0 

0.0078 

1.11 

10 , 

578 

74.0 

11.7 

0.0087 

1.06 

20 

1010 

46.0 

9.5 

0.0101 

0.978 

40 

1595 

22.5 

8.1 

0.0106 

0.909 






Miami Silt 

Loam, 13.0% 

M.C. , 0.2% Salt 


f 

x*f 

R 

e 

O 

Log e 

5 

67 

46.5 

20.9 

0.0207 

1.32 

10 

169 

40.0 

16.2 

0.0220 

1.20 

20 

396 

30.1 

13.5 

0.0237 

1.13 

40 

830 

17. 0 

12.3 

0.0226 

1.09 



Miami silt 

Loam, 15.7% 

M.C. , 0.2% Salt 


f 

x*f 

R 

e 

a 

Log e 

5 

75 

48.0 

21.6 

0.0198 

1.33 

10 

184 

38.7 

18.1 

0.0219 

1.25 

20 

409 

27.9 

15.2 

0.0238 

1.18- 

40 

780 

14.4 

14.2 

0.0232 

1.15 
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Miami Silt 

1 

o 

17. 8% M.C. , 0.2% Salt 


f 

X*f 

R 

e 

a 

Log £ 

5 

148 

57.5 

25.9 

0.0143 

1.41 

10 

307 

48.6 

16.8 

0.0153 

1.22 

20 

677 

35.0 

12.6 

0.0151 

1.10 

40 

1000 

18.2 

11.2 

0.0184 

1.04 




Miami Silt 

Loam/ 

18.2% M.C. / 0.2% Salt 


f 

X*f 

R 

e 

a 

Log E 

5 

110 

45.0 

32. 3 

0.0187 

1.50 

10 

210 

33.0 

25.0 

0.0223 

1.39 

20 

360 

24.0 

17,9 

0.0270 

1.25 

40 

680 

13.1 

15.7 

0.0266 

1.19 




Miami Silt 

Loam/ 

19.1% M.C., 0.2% Salt 


f 

x*f 

R 

e 

a 

Log e 

5 

425 

99.0 

18.1 

0.0061 

1.25 

10 

470 

66.2 

12.7 

0.0104 

1.10 

20 

1100 

30.0 

12.4 

0.0078 

1.09 

40 

1236 

9.55 

12.6 

0.0087 

1.10 
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Miami Silt Loam, 22.0% M.C. , 0.2% Salt 


f 

x*f 

R 

■ E 

a 

Log E 

5 

25 

29.0 

21.0 

0.0349 

1.32 

10 

84 

26.5 

19.7 

0.0355 

1.29 

20 

238 

21.0 

18.2 

0.0364 

1.26 

40 

570 

12.4 

17.0 

0.0322 

1.23 




Miami Silt 

Loam, 10% M.C. , 

0.4% Salt 


f 

x*f 

R 

e 

a 

Log e 

5 

86 

54.0 

19.5 

0.0175 

1.29 

10 

189 

42.4 

15.8 

0,0204 

1.20 

20 

416 

31.1 

13.2 

0.0227 

1.12 

40 

870 

37.2 

12.0 

0,0242 

1.08 




Miami Silt 

Loam, 12.3% M.C. , 

0.4% Salt 


f 

X*f 

R 

e 

0 

Log e 

5 

74 

44.0 

25.1 

0.0213 

1.40 

10 

162 

35.8 

19.0 

0.0240 

1.28 

20 

361 

28.2 

14.3 

0.0256 

1.15 

40 

770 

16.0 

13.1 

0.0243 

1.11 
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Miami Silt 

Loam, 

15.9% H.C. , 0.4% Salt 


f 

x*f 

R 

c 

a 

Log e 

5 

189 

45.2 

40.2 

0.0136 

1.60 

10 

228 

30.5 

28.7 

0.0217 

1.45 

20 

372 

21.0 

21.2 

0.0268 

1.32 

40 

610 

12.0 

17.2 

0.0295 

1.23 




Miami Silt 

Loam, 18.9% 

M.C. , 0.4% Salt 


f 

x*f 

R 

e 

a 

Log e 

5 

53 

33.0 

32.5 

0.0286 

1.51 

10 

135 

25.7 

29.3 

0.0315 

1.46 

20 

250 

18.0 

23.3 

0.0376 

1.36 

40 

420 

10.9 

19.4 

0.0434 

1.28 
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Properties of Some Soils in the 250 - 450 MHz Range 

These data were collected using the General Radio 1602 
Admittance Meter. The procedures used to prepare samples emd 
collect data are presented in Chapter 5, Peurt I. Descriptions 
of the soils used in the investigation are summarized in Ap- 
pendix II. The computer program used to convert the data in- 
to permittivities and conductivities is listed in Appendix IV. 
The symbols used in the tables are as follows : 
f = frequency in MHz 
B = susceptance in mhos 
G = conductance in mhos 
e = relative permittivity 
Log e logarithm to base 10 of e 
a “ conductivity in mhos/meter 
M.C. » moisture content 



* = denotes 

samples prepared with 

salt 




Miami (Lab Prepared} , 

9* M.C. 


f 

B 

G 

c 

Log c 

0 

250 

0.0139 

0.0008 

4.97 

0.697 

0.008 

300 

0.0176 

0.0012r 

4.90 

0.690 

0.010 

350 

0.0219 

0.0011 

4.77 

0.679 

0.012 

450 

0.0350 

0.0035 

4.?3 

0.675 

0.014 


I 


igo 


Miami (Lab Prepared), 10,2% H.C . 


f 

B 

G 

e 

Log c 

a 

250 

0.0148 

0.0011 

5.60 

0.748 

0.010 

300 

0.0189 

0.0015 

5.52 

0.742 

0.012 

350 

0.0238 

0.0022 

5.45 

0.736 

0.014 

450 

0.0392 

0.0046 

5.37 

0.730 

0.016 



Miami (Lab Prepared) ,12.6% 

M.C. 


f 

B 

G 

c 

Log c 

a 

250 

0.0161 

0. 0017 

6.45 

0. 810 

0 . 015 

300 

0.0207 

0.0023 

6.37 

0.804 

0.018 

350 

0.0263 

0.0032 

6.25 

0.796 

0.019 

450 

0.0451 

0.0072 

6.16 

0,790 

0.021 



Miami (Lab Prepared), 14.0% 

M.C. 


f 

B 

G 

e 

Log e 

0 

250 

0.0182 

0.0015 

7.76 

0.890 

0.013 

300 

0.0237 

0.0023 

7.65 

0.884 

0.016 

350 

0.0313 

0.0033 

7.66 

0.885 

0.017 

450 

0.583 

0.0095 

7.50 

0.875 

0.020 
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Miami (Lab Prepared)^ 15.0% M.C . 


f 

B 

G 

e 

Log e 

0 

250 

0.0209 

0.0014 

9.34 

0.970 

0.011 

300 

0.0277 

0.0021 

9.20 

0.964 

0.013 

350 

0.0369 

0.0041 . 

9.04 

0.956 

0.018 

450 

0.792 

0.0117 

8.92 

0.950 

0.016 



Miami 

(Lab Prepared), 17.6% 

M.C. 


f 

B 

G 

c 

Log e 

a 

250 

0.0239 

0.0023 

11.0 

1.04 

0.017 

300 

0.0322 

0.0035 

10.8 

1.03 

0,019 

350 

0.0446 

0.0058 

10.6 

1.02 

0.021 

450 

0.114 

0.0304 

10.4 

1.02 

0.02 3 



Miami (Lab Prepared), 18.8% M.C. 


f 

B 

G 

e 

Log e 

a 

250 

0.0270 

0.0028 

12.6 

1.10 

0.019 

300 

0.0372 

0.0044 

12.3 

1.09 

0,021 

350 

0.0532 

0.0078 

12.1 

1.08 

0.023 

450 

0.169 

0.0686 

11.9 

1.07 

0.025 
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Miaiai (Lab Prepared) , 21«0% M.C . 


f 

B 

G 

E 

Log £ 

a 

250 

0.0341 

0.0043 

15.8 

1.20 

0.025 

300 

0.0499 

0.0075 

15.5 

1.19 

0.027 

350 

0.0802 

0.0161 

15.3 

1.18 

0.028 

450 

-0.178 

0.485 

15.2 

1.18 

0.030 



Miami 

(Lab Prepared), 22.0% 

M.C. 


f 

B 

G 

c 

Log c 

0 

250 

0.0481 

0.0053 

20.9 

1.32 

0.023 

300 

0.0798 

0.0120 

20.6 

1.31 

0.025 

350 

0.173 

0.0492 

20.1 

1.30 

0.026 

450 

-0.148 

0.0265 

19.8 

1.29 

0.028 






Miami (Lab Prepared), 24.2% 

M.C. 


f 

3 

G 

£ 

Log E 

a 

250 

0.0579 

0.0069 

23,7 

1.37 

0.025 

300 

0.107 

0.0192 

23.4 

1.37 

0,027 

350 

0.315 

0.184 

23.2 

1.36 

0.029 

450 

- 0.105 

0.118 

22.9 

1.36 

0.031 
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Miami (Field, 

3" t) 8 

.6% H.C. 


f 

B G 

e 

Log e 

a 

250 

0.0134 0.0010 

4.59 

0.662 

0.010 

300 

0.0169 0.0021 

4.55 

0.658 

0.018 

350 

0.0211 0.0031 

4.51 

0.654 

0.022 

450 

0.0328 0.0076 

4.48 

0.651 

0.032 

t refers to depth below surface of 
caken. 

earth at Which 

sample was 



Miami 

. (Field, 

3"), 11.5% 

M.C. 


f 

B 

G 

c 

Log e 

a 

250 

0.0153 

0.0022 

5.90 

0.771 

0.020 

300 

0.0195 

0.0037 

5.84 

0,766 

0.029 

350 

0.0246 

0.0053 

5.78 

0.762 

0.033 

450 

0.0407 

0.0120 

5.75 

0,760 

0.039 


Miami 

. (Field, 

3"), 17.0% 

M.C. 


f 

B 

G 

. c 

Log c 

a 

250 

0.0258 

0.0043 

12.0 

1.08 

0.030 

300 

0.0356 

0.0074 

11.9 

1.07 

0.037 

350 

0.0504 

0.0136 

11.8 

1.07 

0.042 

450 

0.118 

0.102 

11.7 

1.07 

0.051 
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Miami (FieldyB") , 19.4% M.C. 


f 

B 

G 

E 

Log c 

o 

250 

0.0283 

0.0092 

13.4 

1.12 

0.060 

300 

0.0394 

0.0144 

13.3 

1.12 

0.064 

350 

0.0564 

0.0259 

13.2 

1.12, 

0.067 

450 

0.0641 

0.180 

13.2 

1.12 

0.070 




Micuni (Field, 

3") , 22.6% 

M.C. 


f 

U 

G 

c 

Log e 

a 

250 

0.0423 

0.0172 

19.5 

1.29 

0.081 

300 

0.0640 

0.0337 

19.3 

1.28 

0.085 

350 

0.0919 

0.0911 

19.0 

1.28 

0.090 

450 

-0.110 

0.0704 

19.0 

• 1.27 

0.096 




Miami (Field, 

3"), 23.2% 

M.C. 


f 

B 

G 

e 

Log, e 

a 

250 

0.0557 

' 0.0215 

23.7 

1.37 

0.078 

300 

0.0928 

0.0566 

23.6 

1.37' 

0.084 

350 

0.0607 

0.222 

23.4 

1.37 

0.094 

450 

- 0.0873 

0.0268 

23.4 

1.37 

0.102 
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Miami (Field, 3")^ 24.Qft M. C. 


f 

B 

G 

c 

Log e 

a 

250 

0.0530 

0.0338 

24.0 

1.38 

0.122 

300 

0.0726 

0.0795 

23.8- 

1.37 

0.130 

350 

0.0020 

0.177 

23.6. 

1.37 

0.134 

450 

-0.0785 

0.0301 

23.4 

1.37 : 

0. 140 




Miami (Field, 

3") , 24.0% M.C. 



f 

B 

G 

e 

Log e 

q 

250 

0.0484 

0.0348 

22.9 

1.36 

0.135 

300 

0.0615 

0.0755 

22.7 

1.35 

0.146 

350 

0.0084 

0,150 

22.5 

1.35 

0.152 

450 

— 0.0766 

0.0357 

22.5 

1.35 

0.158 



Miami (Fields 18**), 8.B% M.C. 


f 

B 

G 

E 

Log e 

a 

250 

0.0135 

0.0009 

4.68 

0.670 

0.009 

300 

0.0171 

0.0016 

4.63 

0.666 

0.013 

350 

0.0213 

0.0034 

4.59 

0.662 

0.024 

450 

0.0335 

0.0073 

4.58 

0.661 

0.030 


I 



Miami 

(Field, 18"), 

10.0% M.C. 

' 


£ 

‘ B 

G 

e 

Log e 

a 

250 

0.0141 

0.0022 

5.13 

0.710 

0.021 

300 

0.0179 

0.0029 

5.08 

0.706 

0.024 

350 

0.0224 

0.0049 

5.05 

0.703 

0.033 

450 

0.0356 

0.0102 

5.01 

0.699 

0.039 



Miami 

(Field, 18"), 

15.0% M.C. 



f 

£ 

G 

E 

Log E 

a 

250 

0.0174 

0.0067 

7.43 

0.871 

0.057 

300 

0.0224 

0.0090 

7.38 

0.868 

0.063 

350 

0.0264 

0.0134 

7.34 

0.866 

0.072 

450 

0.0444 

0.0339 

7.33 

0.865 

0.084 



w 


Miami (Field, 18"), 16.0% M.C. 


f 

. B 

G 

e 

Log £ 

a 

250 

0.0217 

0.0074 

9.98 

0.999 

0.079 

300 

0.0288 

0.0108 

9.93 

0.997 

0.088 

350 

0.0379 

0.0175 

9.84 

0.993 

0.100 

450 

0.0617 

0.0627 

9.77 

0.990 

0.109 



Miami 

(Field, 

18") , 17.8% M.C. 



f 

B 

G 

- £ 

Log £ 

a 

250 

0.0244 

0.0088 

11.4 

1.06 

0.063 

300 

0.0327 

0.0142 

11.4 

1.05 

0.074 

350 

0.0446 

0.0239 

11.3 

1.05 

0.081 

450 

0.0573 

O.iOO 

11.3 

1.05 

0.092 



Miami 

(Field, 

18") , 18.6% M.C. 



f 

B 

G 

E 

Log £ 

0 

250 

0.0301 

0.0127 

14.4 

1.16 

0.079 

300 

0.0413 

0.0213 

14.2 

1.15 

0.088 

350 

0.0554 

0.0425 

14.2 

1.15 

0.100 

450 

-0.0199 

0.146 

14.1 

1.15 

0. 109 


Miami (Field, 18"), 20.4% M.C. 


f 

B 

G 

e 

Log e 

a 

250 

0.0358 

0.0179 

17.1 

1.23 

0.096 

300 

0.0499 

0.0322 

16.9 

1.23 

0.105 

350 

0.0606 

0.0721 

16.9 

1.22 

0.120 

450 

-0.0798 

0.101 

16.7 

1.22 

0.125 



Micuoni (Field, 

18") , 22.0% M.C. 



f 

B G 

E 

Log c 

a 

250 

0.0385 0.0231 

18.6 

1.27 

0.115 

300 

0.0530 0.0429 

18.4 

1.26 

0.123 

350 

0.0590 0.0921 

18. 3 

1.26 

0.128 

450 

-0.0849 0.0759 

18.2 

1.25 

0.135 


Miami (Field, 

18") , 23.0% M.C. 




f 

B 

G 

e 

Log £ 

0 

250 

0.0433 

0.0250 

20.4 

1.31 

0.112 

300 

0.0621 

0.0510 

20.4 

1.31 

0.119 

350 

0.0591 

0.126 

20.4 

1.31 

0.124 

450 

-0.0887 

0.0507 

20.4 

1.31 

0.130 



1S9 


Miami (Field, 18"), 23.18% M.C. 


f 

B 

G 

e 

Log c 

a 

250 

0.0457 

0.0299 

21.6 

li33 

0,125 

300 

0.0629 

0.0623 

21.5 

1.33 

0,132 

350 

0.0344 

0.141 

21.4 

1.33 

0.140 

450 

-0.0813 

0.0426 

21.3 

1.33 

0.147 
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Miami (Field, 36”)/ 9»0% M,C. 


f 

B 

G 

e 

Log e 

a 

250 

0.0135 

0.0012 

4.67 

0.670 

0.012 

300 

0.0171 

0.0023 

4.63 

0.666 

0.019 

350 

0.0224 

0.0036 

4.59 

0.662 

0.025 

450 

0.0333 

0.0077 

4.56 

0.659 

0.032 


Miami 

(Field, 

36") , 11.1% 

M.C. 


f 

B 

G 

£ 

Log e 

a 

250 

0.0139 . 

0.0020 

5.01 

0.710 

0.019 

300 

0.0178 

0.0033 

5.01 

0.710 

0.027 

350 

0.0224 

0.0048 

5.01 

0.710 

0.032 

450 

0.0357 

0.0100 

5.01 

0.710 

0.038 


Miami 

(Field, 

36") , 14.0% 

M.C. 



f 

B 

G 

E 

Log e 

a 

250 

0.0171 

0.0034 

7.08 

0.850 

0.030 

300 

0.0219 

0.0053 

6.99 

0.845 

0.038 

350 

0.0282 

0.0079 

6.97 

0.843 

0.044 

450 

0.0483 

0.0201 

6.90 

0.839 

0.051 
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Miami 

. (Field, 36"), 

18.6% M.C. 



f 

B 

G 

e 

Log e 

a 

250 

0.0251 

0.0098 

11. 8 

1.07 

0.069 

300 

0.0337 

0.0146 

11.7 

1.07. 

0.074 

350 

0.0457 

0.0235 

11.6 

1.06 

0.077 

450 

0.0712 

0.103 

11.4 

1.06 

0.080 


f 

Miami 

B 

(Field, 

G 

36"), 

21.0% M.C. 
e 

Log e 

a 

250 

0.0392 

0.0138 


18.2 

1.26. 

0.070 

300 

0.0580 

0.0270 


17.9 

1.25 

0.078 

350 

0.0881 

0.0690 


17.9 

1.25 

0.083 

450 

-0.118 

0.0967 


17.8 

1.25 

0.088 

f 

Miami 
. B 

(Field, 

G 

36"), 

22.2% M.C. 

e 

Log c 

a 

250 

0.0479 

0.0313 


22.3 

1.35- 

0.125 

300 

0.0661 

0.0658 


22.1 

1.34 

0.130 

350 

0.0339 

0.150 


21.8 

1.34 

0.135 

450 

-0.0829 

0.0408 


21.6 

1.33 

0.140 
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Chelsea (Lab Prepared) , 

2.3% M.C. 

f 

B 

G 

e 

Log e 

250 

0.0128 

0.0010 

4.19 

0.622 

300 

0.0161 

0.0014 

4.13 

0.616 

350 

0.0200 

0.0018 

4.10 

0.613 

450 

0.0310 

0.0031 

4.05 

0.607 


Chelsea (Lab Prepared), 4.04 M.C. 


f 

B 

G 

e 

Log £ 

250 

0.0136 

0.0009 

4.77 

0.679 

300 

0.0173 

0.0013 

4.75 

0.677 

350 

0.0217 

0.0019 

4.72 

0.674 

450 

0.0347 

0.0035 

4.68 

0.671 



Chelsea (Lab 

Prepared ) , 

5.4% M.C. 

f 

B 

G 

e 

Log e 

250 

0.0142 

0.0016 

5.19 

0.715 

300 

0.0182 

0.0020 

5.16 

0.713 

350 

0. 0229 

0.0029 

5.15 

0.712 

450 

0.0372 

0.0054 

5.11 

0.708 


0.010 

0.012 

0.013 

0.014 


a 

0.009 

0.0105 

0.013 

0.014 


a 

0.015 

0.0165 

0.019 

0.020 
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Chelsea (Lab Prepared) ,7.0% M.C . 


f 

B 

G 

e ; 

Log e 

250 

0.0159 

0.0021 

6.34 

0.802 

300 

0.0204 

0.0030 

6.25 

0.796 

350 

0.0261 

0.0045 

6.22 

0.794 

450 

0.0447 

0.0106 

6.19 

0.792 



Chelsea (Lab Prepared) , 

8.6% M.C. 

f 

D 

G 

e 

Log E 

250 

0.0164 

0.0025 

6.61 

0.820 

300 

0.0212 

0.0033 

6.61 

0. 820 

350 

0.0274 

0.0045 

6.61 

0.820 

450 

0.0491 

0.0107 

6,66 

0.824 



Chelsea (Lab Prepared) , 

9.1% M.C. 

f 

B 

G 

£ 

Log c 

250 

0.0177 

0.0039 

7.50 

0.875 

300 

0.0231 

0.0053 

7.46 

0.873 

350 

0.0300 

0.0074 

7.43 

0.871 

4S0 

0.0542 

0.0186 

7.41 

0.870 




a 

.0,019 

0.023 

0.027 

0.031 


a 

0.022 

0.024 

0.026 

0.028 


a 

0.034 

0.037 

0.039 

0.041 
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Chelsea (Lab 

Prepared) . 

10.1% M.C, 


f 

B 

G 

c 

Log c 

o 

250 

0.0185 

0.0045 

7.99 

0.903 

0.038 

300 

0.0242 

0.0059 

7.94 

0.900 

0.040 

350 

0.0318 

0.0084 

7.94 

0.900 

0.042 

450 

0.0603 

0.0226 

7.96 

0.901 

0.043 
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Crider (Lab Prepared) , 14.5% M.C. 


f 

B 

G 

c 

Log e 

a 

250 

0.0130 

0.0031 

4.37 

0.640 

0.030 

300 

0.0163 

0.0035 

4.25 

0.628 

0.030 

350 

0.0197 

0.0069 

4.17 

0.620 

0.050 

450 

0.0290 

0.0130 

4.07 

0.610 

0.060 



Crider 

(Lcdi Prepared) , 16.0% 

M.C. 


f 

B 

G 

e 

Log e 

0 

250 

0.0133 

0.0062 

4.69 

0.671 

0.060 

300 

0.0164 

0.0083 

4.56 

0.659 

0.070 

350 

0.0201 

0.0099 

4.48 

0.651 

0.070 

450 

0.0287 

0.0179 

4.36 

0.640 

0.080 



Crider 

(Lab Prepared), 17.0% 

M.C. 


f 

B 

G 

e 

Log e 

0 

250 

0.0128 

0.0125 

4.76 

0.678 

0.120 

300 

0.0152 

0.0155 

4.68 

0.670 

0.130 

350 

0.0171 

0.0194 

4.52 

0.655 

0.140 

450 

0.0195 

0.0299 

4.41 

0.644 

0.150 
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Crider 

(Lab Prepared) r 18.6% 

M.C. 


f 

B 

G 

e 

Log E 

a 

250 

0.0141 

0.0107 

5.48 

0.739 

0.100 

300 

0.0169 

0.0137 

5.31 

0.725 

0.110 

350 

0.0203 

0.0164 

5.20 

0.716 

0.110 

450 

0.0262 

0.0284 

5.07 

0.705 

0.120 



Crider 

(Lab Prepared), 21.6% 

M.C. 


f 

B 

G 

e , 

Log E 

a 

250 

0.0146 

0.0120 

5.95 

0.775 

O.llO 

300 

0.0176 

0.0154 

5. 83 

0.766 

0.120 

350 

0.0204 

0.0202 

5.73 

0.758 

0.130 

450 

0.0241 

0.0346 

5.62 

0.750 

0.140 



♦Crider 

' (Lab Prepared) , 

22.5% M.C. 


f 

B 

G 

E 

Log C 

a 

250 

0.0157 

0.0174 

7.24 

0.860 

0.150 

300 

0.0182 

0.0221 

7.05 

0.848 

0.160 

350 

0.0207 

0.0272 

6.93 

0.841 

0.160 

450 

0.0185 

0.0455 

6.76 

0.831 

0.170 
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Crider (Lab Prepared) 


25.0% M.C. 


f 


G 

e 

Log £ 

a 

250 

0.0185 

0.0134 

8.55 

0.932 

0.110 

300 

0.0227 

0.0181 

8.29 

0.919 

0.120 

350 

0.0267 

0.0255 

8.14 

0.911 

0.130 

450 

0.0284 

0.0516 

7.94 

0.899 

0. 140 




Crider (Lab Prepared), 26.5% 

M.C. 


f 


G 

e 

Log e 

o 

250 

0.0192 

0.0143 

8.99 

0.954 

0.115 

300 

0.0237 

0.0187 

8.75 

0.942 

0.120 

350 

0.0280 

0.0267 

8.61 

0.935 

0.130 

450 

0.0286 

0.0561 

8.41 

0.925 

0.140 




*Crider (Lab Prepared) , 

27. 4% M.C. 


f 

B 

G 

e 

Log e 

a 

250 

0.0196 

0.0180 

9.59 

0.982 

0.140 

300 

0.0232 

0.0242 

9. 31 

0.969 

0.150 

350 

0.0257 

0.0336 

9.14 

0.961 

0.160 

450 

0.0120 

0.0627 

9.02 

0.955 

0.180 
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Crider (Lab Prepared), 29.2% M.C, 


f 

B 

G 

E 

Log c 

0 

250 

0.0219 

0.0207 

11.1 

1.04 

0.150 

300 

0.0258 

0.0287 

10. 8 

1.03 

0. 160 

350 

0.0276 

0.0413 

10.7 

1.03 

0.170 

450 

0.0085 

0.0759 

10.4 

1.02 




*Crider (Lcib Prepared) ^ 30.6% M.C. 


f 

B G 

e Log E 

a 

250 

0.0234 0.0216 

12.1 1.08 

0.150 

300 

0.0285 0.0313 

12.0 1.08, 

0.160 

350 

0.0303 0.0476 

12.0 1.08 

Q.170 

450 

.0.0092 0.0843 

12.0 1.08, 

0.180 


Crider (Lab Prepared) 32.0% M.C. 



f 

B 

G 

E 

Log E 

0 

250 

0.0298 

0.0216 

15.0 

1.17 

0.130 

300 

0.0374 

0.0341 

14.6 

1.16 

0.140 

350 

0.0442 

0.0555 

14.4 

1.15 

0.140 

450 

^-0.0305 

0.106 

14.1 

1.15 . 

0.150 
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Crider (Lab Prepared), 34.4% M.C. 


f 

B 

G 

e 

Log £ 

a 

250 

0.0318 

0.0291 

16.8 

1.22 

0. 160 

300 

0.0373 

0.0463 

16.3 

1.21 

0. 170 

350 

0.0307 

0.0756 

16.0 

1.20 

0.180 

450 

-0.0512 

0.0822 

15.6 

1.19 

0.180 



Crider 

' (Lab Prepared), 36.4% 

M.C. 


f 

B 

G 

£ 

Log e 

0 

250 

0.0354 

0.0345 

18.9 

1.27 

0.170 

300 

0.0398 

0.0576 

18.4 

1.26 

0.180 

350 

0.0211 

0.0938 

18.1 

1.25 

0.190 

450 

-0.0619 

0.0614 

17.7 

1.25 

0.200 



Crider (Lad> Prepared), 36.4% 

M.C. 


f 

B 

G 

e 

Log c 

a 

250 - 

0.0421 

0.0376 

21.4, 

1.33 

0.160 

300 

0.0476 

0.0688 

20.8 

1.31 

0.170 

350 

0.0193 

0.119 

20.4 

1.31 

0.170 

450 

-0.0697 

0.0491 

19.9 

1.30 

0.190 


Crider (Lab Prepared ) , 40.0% M.C . 


f 

B 

G 

e 

Log c 

a 

250 

0.0513 

0.0465 

25.0 

1.39 

0.160 

300 

0.0549 

0.0956 

24.3 

1.38 

0.165 

350 

-0.0227 

0.149, 

23.8 

1.37 

0.166 

450 

-0.0726 

0.0313 

23.4 

1.37 

0.170 



201 



Crosby 

(Field, 

18"), 14.0% 

M.C. 


f 

B 

G 

e 

Log e 

a 

250 

0.0139 

0.0074 

5,13 

0.710 

0.070 

300 

0,0171 

0.0098 

5.04 

0.702 

0.080 

350 

0.0209 

0.0122 

4.98 

0.697 

0.083 

450 

0.0299 

0.0221 

4.91 

0.691 

0.091 



Crosby 

(Field, 18") 

, 16.0% 

M.C. 


f 

B 

G 

£ 

Log e 

a 

250 

0.0173 

0.0063 

7.33 

0.865 

0.055 

300 

0.0221 

0.0092 

7.23 

0.859 

0.065 

350 

0.0279 

0.0132 

7.18 

0.856 

0.072 

450 

0.0445 

0.0299 

7.08 

0.850 

0.077 



Crosby 

(Field, 

18"), 19.8% 

M.C, 


f 

B 

G 

e 

Log e 

a 

250 

0.0198 

0.0113 

9.12 

0.960 

0.090 

300 

0.0249 

0.0158 

8.91 

0.950 

0-100 

350 

0.0312 

0.0223 

8.81 

0.945 

0.104 

450 

0.0406 

0.0553 

8.71 

0.940 

0.113 



2Q2 



Crosby 

(Field, 

18"), 21.2% 

M.C. 


f 

B 

G 

e 

Log z 

o 

250 

0.0236 

0.0146 

11.4' 

1.06 

0. 105 

300 

0.0305 

0.0214 

11.3 

1.05 

0.113 

350 

0.0377 

0.0339 

11,2 

1.05 

0.122 

450 

0.0255 

0.0922 

11.1 

1.04 

0.129 



Crosby 

(Field, 

18") , 24% M.C. 



f 

B 

G 

e 

Log £ 

0 

250 

0.0323 

0.0231 

16.2 

1.21 

0.130 

300 

0.0412 

0.0380 

15.9 

1.20 

0.140 

350 

0.0464 

0.0664 

15.7 

1.19 

0.143 

450 

-0. 0524 

0.0998 

15.4 

1.18 

0.151 



Crosby 

(Field, 

18"), 25.8% 

M.C. 


f 

B 

G 

e 

Log e 

or 

250 

0,0370 

0.0294 

18.7 

1,27 

0.145 

300 

0.0474 

0.0514 

18.6 

1.27 

0.151 

350 

0.0404 

0.0976 

18.4 

1.26 

0.155 

450 

-0.0751 

0.0671 

18.3 ■ 

1,26 

0.160 
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Crosby (Fields 18”), 26.8% M.C. 


f 

B 

G 

■■ e' 

Log e 

a 

250 

0.0377 

0.0327 

19.4 

1.28 

0.156 

300 

0.0463 

0.0573 

19.2 

1.28 

0.162 

350 

0.0312 

0.104 

19.0 

1.28 

0.166 

450 

-0.0717 

0.0595 

18.8 

1.27 

0.175 
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APJ>ENDIX II 

DESCRIPTION OF SOME SOILS USED IN THIS INVESTIGATION 

Chelsea 

Chelsea is a moderately dark colored soil which is deep 
and well-drained. The surface of the soil (where samples for 
this investigation were collected) is dark brown, while the 
subsoil is yellowish red to reddish brown. The soil is more 
than 80% sand and contains very little organic matter. The 
wilting point moisture for the semples used in this investi- 
gation was 2.2%. 

Crider 

Where it occurs, Crider is a deep well-drained, gently 
sloping soil. Its surface layer consists of dark brown silt 
loam while the subsurface (where sounples for this investi- 
gation were obtained) is silty clay or clay which is red. 
Crider contains between 27% and 40% clay. The wilting point 
of the samples used in this investigation was 18.5%. 

Crosby 

Crosby developed under hardwood forests and occurs on 
nearly level ground. It is somewhat poorly drained and dark 
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in color. On the surface there is a layer of silt loam and 
below that is brown silt loam (where seunples for this in- 
vestigation were taken) . The subsoil consists of thick yel- 
lowish brown silty clay loam. The wilting point for samples 
used in this investigation was 14%. 

Miami 

Miami occurs in upland tills and moraines on sloping 
ground. It is well-drained, deep, and moderately dark in 
color. The soil is predominately silt, more than 66%. In 
the subsoil, however, both sand and clay content increase 
somewhat. The moisture wilting point for the samples used 
in this investigation was 6.8%. 
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APPENDIX III 

I 

ANALYSIS OF SAMPLE HOLDER-SOIL EQUIVALENT CIRCUIT 


The circuit is shown in Figure A-1. The impedance of 
this circuit is 


R, 


A-1 Z 


jo)C, 


jwC, 


12 


-= + Rj^ -= — + R2 


^1 ^ ^2 




'2 *'2 


A-2 Z,, = 


Rj^^ - jwCj^Rj^* R2 - 3“C2R2* 


1 + u*Cj^®Rj^* 1 + u*C2*R2^ 
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R| = Electrode resistance 
C| = Electrode capacitance 
R2 = Sample resistance 
C2 = Sample copacitance 
R > Measured resistance 
C ~ Measured capacitance 

Figure A-1. Sample Soil Holder Circuit. 
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A- 3 


'12 


(1 + a)*Cj^®Rj^*} (1 + u*C2*R2M 

j(o[Cj^Rj^Ml + w*C 2 *R 2 *> + C2R2M1 + u^Cj^^Rj^*)] 
(1 + u*Cj^*Rj^M (1 + w*C2*R2^^ 


A_4 Z - --- 


1 + (J^R^C® 


A- 5 


1 + w*R^C^ 


Rj^d + u)^C2*R2^) + R 2U + 

(1 + w^Cj^^Rj^^) (1 + 64^C2^R2^) 


Equating the imaginary part of Z to the imaginary part of Z^2 


A- 6 


CR^ 


1 + u*R*C* 


CiRi^ (1 + U)^C2^R2M + C2R2^(1 + 

(1 + (1 + uj®C2*R2*> 


From A- 5 


A-7 1 + RC = 


R(1 + (1 + u*C2^R2*^ 

R^(l + (u^C2*R2M + R2(1 + 


Substituting into A-6 


Cj^Rj^MI + (u*C2*R2^) + C 2 R 2 MI + 

Rj^d + u)^C2*R2^) + R 2 U + u)*C^^R^M 


A-8 RC = 
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Substituting into A- 5 


A-9 R 




Rj^d + a)*C2*R2^) ''' ^2 U)*Cj^*Rj^M 


Rj^d + 0 )*C 2 *R 2 *) + OJ^Cj^^Rj^*) 

(1 + !u®C2*R2*) 


Simplifying 


Rj^MI + u*R 2“C2^) + 2Rj^R2d + u ^ Cj^Rj^C2R2) + R2^ (1 + 


(1 + + R2(1 + (J^C^^Rj^^) 


A-ll c = |£ = 


CiRi* (1 + u)^C 2^R2^) + C2R2^(1 + t*>* Rj^^ ) 

Rj^MI + u*R 2*C2*) + 2R^R2(1 + w*Cj^Rj^C2R2) + R2^ U +w * Cj^ ^ R^^ ^ } 


+(cJRj^^+C^R2^):+Cj^C2(Ri+ R^/ [Cj^C2 {€^^ 4 - 02 ) "Ri^R2' 

C]^ C 2 uj* + (Rj^ + R2 )V[(Cj^ + C2)*R]_®R2*] 
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' '*' C2)*Ri^R2*^ 

Rj^Cj^* + R 2 ^ 2 ^ “* * ^*^1 ^*^ 1*^2 ^^ 2 ^ 2 ^ ''■ 


The final results are 


A-1 4 C = 


'^l''2 


^ CiRi^ 

CiC2<Ci + C2 )Rj^*R2® 


C, + C, 

^ 2 »" + 


(R^ + R2)* 


<Ci + C2 )^Ri*R2* 


A-15 R = 


RiR2(Ci ->• C^)^ 

RlCi* + R2C2* 


U)* + 


+ 


(Rj^ + R2)* 


(Cj^ + C 2 ) *Ri*R2* 


®i * ■'a 


Va'^a'^a' * Vi’’ 




